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STANISLAW LORIA (1883—1958) 
I: 


Professor Loria died on the 8th of August at the age of 75 in London. Before I say 
:more about his life and work I should like to give a few personal reminiscences. 
I met Professor Loria when I was a student in Kraków. He was then a young docent, 
:and I was just starting my studies. He lectured on theoretical mechanics, of which 

I knew almost nothing. Although I did not attend all the lectures, since they were 
too difficult for me, still the few lectures I attended remain vividly in my mind, because 
of the beauty of the subject and of his presentation. During the second year of my 
: studies I met Docent Loria at the physic’s lab. And I well remember his helpful 
attitude to young, eager students — an attitude full of grace and charm. . 

Later, when Professor Smoluchowski died, Stanislaw Loria went to Lwöw to replace 
Professor Zakrzewski, who moved into Smoluchowski's chair. For years I lost contact 
with Professor Loria. As a high school teacher in the provinces and then in Warsaw, 

. I had hardly any opportunity to see him. We met again at the Physics Meeting in 
Lwöw, the first which I attended and at which I presented my first paper. I was very 
much touched by the fact that Professor Loria remembered me and by his attitude 
towards me. He became interested in my work and he asked me to send him my 
manuscript. The next time I saw Professor Loria was some years later in Berlin in 
1928 during a course in theoretical physics at which Schròdinger and Laue lectured. 
I think that it was there that our fiendship started. It was Professor Loria who later” 
offered me an assistantship at Lw6w University and it was he who took me out of the 
high school and started my academic career. Indeed, I do not owe to anyone else as 
much as I owe to Professor Loria. I have written this introduction to lead up to this 
statement which I should like to record now that he is not among the living. 

There is an additional reason for my writing this. For it shows his attitude towards 
younger people, his care for their well-being and scientific development. These are 
the characteristic features of a good teacher. And a good teacher he was! 


IL. 


| 
| 
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Stanislaw Loria was born in 1883 in Warsaw. He went to elementary school, high 
school and university in Kraköw. After his Ph. D. he went for further studies to 
Wroclaw, Gòttingen, Berlin. In the year 1910 he became a docent (reader) at the 
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Jagiellonian University in Kraków. Three years later he visited England and worked 
at the University of Manchester. When Professor Zakrzewski, who had the chair of 
theoretical physics in Lw6w, went to Krakéw, Stanislaw Loria became a professor in 
Lwöw to replace him. The years 1923—1926 he spent in Pasadena as a research associate 
in the California Institute of Technology. When he returned to Poland in 1927, he 
occupied the chair of experimental physics at Lwöw University. There he remained 
until 1941, that is until the Germans occupied the town. He had to hide during the 
subsequent four years of occupation. When Poland was freed, he became the vice- 
-delegate of the University of Education for Lower Silesia, then a professor and vice- 
-rector at the University of Wroclaw. In the years 1947—1948 he received a UNESCO 
fellowship to England. In 1951 he moved from Wroclaw to Poznan as a professor 
of experimental physics and as a scientific worker in the Institute of Physics of the 
Polish Academy of Sciences. He died on the 8th of August 1958 in a London 
hospital, where he had undergone an operation. 

The death of Stanislaw Loria is a great loss to Polish physics. He was one of our 
oldest professors, one who knew well almost all the great physicists of the first half of 
our century. He knew Lorentz and Ehrenfest, studied with Born and Ladenburg in 
Breslau. (Usually Professor Born’s first words to me, whenever we met were: “How 
is Loria?"). He worked in Manchester under Lord Rutherford and for the years 1923— 
1926 he was invited by professor Millikan to Pasadena. He was well acquainted with 
the old centers of research. It was a delight to hear him talking about his great teacher 
Witkowski, about Smoluchowski and Natanson, professors of physics in the Univer- 
sity of Krak6w, about Madame Sktodowska Curie, about all the great physicists whom 
he knew personally and whom the next generation knew only from literature. 


II. 


Professor Loria’s scientific work was concentrated around five problems which 
we shall mention briefly: 
First: Papers in the field of experimental psychology: 
1. Untersuchungen über das seitliche Sehen, Bull. Acad. Cracovie, October 1904. 
2. Uber das periphere Sehen, Zeitr. f. Psychol. u. Physiol. d. Sinnesorgane 40, 
160—186, 1906. ' i 
'These papers were important at the time since they proved that one of Wundt's 
philosophical theories concerning “attention” was contradicted by experiment. 
Second: Papers in the field of optics: 
1. Untersuchungen über die Dispersion des Lichtes in Gasen. II. Die Dispersion 
des Acetylens u. Methans, Bull. Acad. Cracovie 1059—1067, 1908. 
2. Uber die Dispersion des Lichtes in leuchtendem Wasserstoff (with R. Laden- 
burg), Nature 79, 7, 1908; Verhandl. d. Deutsch. ck Ges. 10, 878—866, 
1908; Phys. Zeitr. 9, 875—878, 1908. 
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Untersuchungen iiber die Dispersion des Lichtes in Gasen. II. Die Dispersion 


des Athylens u. Athans. Bull. Acad. Cracovie, 195—207, 1909. 


. Über die Disperion des Lichtes in gasformigen Kohlenwasserstoffen, Ann. d. 


Physik, 29, 605—622, 1909. 


. Bestimmung der Dispersion des Lichtes in nicht leuchtendem gesättigtem Na- 


-Dampfe von der Temperatur 385 Gr. Celsius, Bull. Acad. Cracovie 39—60, 
1909; Ann. d. Phys. 30, 240—356, 1909. 


. Untersuchungen über die Dispersion des Lichtes in Gasen. III. Ammoniak 


(with J. Patkowski), Bull. Acad. Cracovie 494—508, 1913. 


. Die Lichtsbrechung in Gasen als physikalisches und chemisches Problem 


(Sammlung Vieweg H. 4.) Fr. Vieweg u. Sohn, Braunschweig 1914. 


. Indirectly Excited Fluorescence Spectra, Phys. Review 26, 573—584, 1925. 
9: 


The Metastable 2p. — State of Mercury Atoms, Proc. Nat. Acad., Washington 
11, 673—679, 1925. 
Zur Frage nach de Abhängigkeit der sensibilisierten Fluorescenze vom Zusatz 


der Gase, Zeitr. f. Physik 38, 672—674, 1926.. 


These papers cover the period of Loria’s work in Wroclaw and Pasadena. The 
Wroclaw papers were the first to show the dispersion of light in some gaseous substan- 
ces. Especially important are the papers written with Ladenburg in which for the first 
time the anormal dispersion in excited hydrogen was discovered. The papers belonging 
to the Pasadena period deal with the problem of fluorescence; they estimated the time 
interval in which the Hg atoms are in a “metastable” state. 


Third: Papers on the magneto-optic Kerr effect: 


L 


4. 


5. 


The Magneto-optic Kerr-effect in Ferromagnetic Compounds and Alloys. I. Proc. 
Amsterdam 12, 835—845, 1910. 


. Uber die optischen Konstanten einiger das magnetooptische Kerr Phanomen 


aufweisenden Substanzen, Bull. Acad. Cracovie (A) 22, 278, 1910. 


. O magneto-optyeznym zjawisku Kerra w ferromagnetycznych związkach i sto- 


pach, Kosmos 35, 621—640, 1910. 

The magneto-optic Kerr-effect in ferromagnetic compounds and alloys. II. 
Proc. Amsterdam 14, 970—983, 1912. 

Der magnetooptische Kerr-Effect bei ferromagnetischen Verbindungén und 
Legierungen, Annalen d. Physik 38, 889—920, 1912. 


These reports on the work done in the laboratory of Professor H. Dubois in Berlin 
were the first in which the Kerr effect in ferromagnetic substances was thoroughly in- 
vestigated. 

Fourth: Papers on radioactivity: 


i 


E 


Ds 


Uber die Verfliichtigung kondensierter Emanationen, Wiener Berichte 124, 


820. 841, 1915. 


Uber die Verdampfung des ThB und ThC, Wien. Beriechte 124, 567—575, 


1915. 
Uber die Verdampfung des RaC, Wien. Berichte 124, 1077—1088, Ds 


5 ; Stanistaw Loria (1883—1958) 


4. Über die Verzwiegungsstelle der Th-Reihe, Phys. Zeitr. 17, 6—9, 1916. 

5. Die Verflüchtigungskurven des Systems ThB + ThC auf Au, Bull. Acad. 
Cracovie, 8/10, 260—264, 1917. 

6. Über die Verflüchtigung aktiver Niederschläge, Bull. Acad. Cracovie 549—556, 
1917. 

These papers indicated the branching out of the radioactive families of Th and Ra, 

Fifth: Work on diffraction of electric waves: 

1. Beugungsversuche mit Elektronen mittlerer Energie, Bull. Acad. Cracovie (A) 

1/2, 15—22, 1937. 

This work forms the first publication on the problem to which Professor Loria 
devoted guite a few years of his life. It deals with diffraction of electrons in the interval 
between 18 and 4.5 thousands volts. Unfortunately, no further publications were 
possible since the apparatus and notes were taken away or destroyed by the Nazis. 

One characteristic of Loria s work would be lacking if Idid not mention his interest 
in theoretical physics and popularisation. Indeed, to give only one example: he was 
one of the first physicists who studied Einstein's paper of 1905 and some years later 
he wrote in his vivid style a book on relativity theory, the first that appeared on this 
subject in Poland. 

Stanisław Loria was a fascinating lecturer, a wonderful teacher and a loyal colleague. 
His scientific work spans the time from the birth of modern physics to the atomic 
age. He was well aware of the scientific problems connected with this rapid development 
and of its political implications. He was wise and decent and had a fine sense of humor. 
Those who knew him will cherish the memory of his charming personality. 


Leopold Infeld 
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DAS VARIATIONSPRINZIP DER RELATIVISTISCHEN 
HYDRODYNAMIK 


Von CZESLAW JANKIEWICZ 


Fakultat fiir Theoretische Physik an der Universitat, Wroclaw; 
Mathematisches Institut der Polnischen Akademie der Wissenschaften, Wroclaw 


(Eingegangen am 23. Februar 1958) 


Auf Grund der Theorie der Raumdeformationen von Schouten und van Kampen 
aus dem Jahre 1953 wurde ein Variationsprinzip fiir die relativistische Hydrodynamik 
aufgestellt und an den Beispielen der idealen und der zähen Flüssigkeit erlaüztert. 


1. Einleitung 


Die Gleichungen der Hydrodynamik bilden ein System partieller Differential- 
gleichungen erster Ordnung. Diese Gleichungen. mit Ausnahme spezieller Fälle 
(z. B. bei inkompressiblen Flüssigkeiten, Potentialstròmungen) lassen sich nicht als 
Euler-Lagrangesche Gleichungen eines entsprechenden Variationsprinzips darstellen. 
Diese Schwierigkeiten kann män beheben durch Einführung von Hilfsfunktionen, 
durch deren Ableitungen sich die Zustandsfunktionen ausdrücken. Zu solchen Var- 
iationsprinzipen gehören die von Bateman (1932), Ito (1953) und Skobelkin (1956). 
Einen anderen Typus stellen die Variationsprinzipen von Fock (1955, $ 47) und Ry- 
barski (1956) dar. Bei diesen variert man den Raum und nicht die Zustandsfunktionen 
(für ein elektrisch geladenes Kontinuum). In der vorliegenden Arbeit wird an Hand 
der Raumdeformationstheorie (Schouten, van Kampen 1933) das Variationsprinzip 
für die relativistische Hydrodynamik formuliert. 


2. Die Raumdeformation 


Wir wollen im Raume X, (s.z.B. Schouten und Struik 1933) die Gruppe der in- 
- finitesimalen Transformationen 

. = 2 + o* (af) de (1) 
betrachten, wo omg (x) ein kontravarianter Vektor und de eine infinitesimale Grósse 
ist (Die durch kleine griechische Buchstaben bezeichneten Indizes nehmen die Werte 
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A; 2, 3, 4 an). Es seien p,(x”) die Komponenten eines geometrischen Obiektes (4 = 
1, 2, ..., N). Den Ausdruck 


Dya = pate) — valv?) ` (2) 


nennen wir das Liesche Differential dieses Obiektes; ‘y ,(x”) bezeichnet hier das durch 
die Transformation (I) deformierte Obiekt a ,(a”). Für geometrische Obiekte der 
Klasse 1, deren Transformationsrelationen nur von den Ableitungen erster Ordnung 
S abhängen, nimmt das Liesche Differential (2) die Form 


Dy, = (W29 Pa — G4,0* — G4, Sa) de (3) 
an (s. z.B. Yano 1955), wo G4,, G4,” für das vorgegebene geometrische Obiekt spezi- 


fische Funktionen sind. Im Folgenden schreiben wir die Lieschen Differentiale fiir 
einige geometrische Obiekte auf. Von diesen Formeln werden wir später Gebrauch 


machen. 
p — Skalar 
Dp = ( pp) de (4) 
o — Skalardichte mit dem Gewicht + 1 
Do = (o^ y,o + eV w") de (5) 
uz — kontravarianter Vektor 
= (op yu" — upa") de (6) 
gas — kontravarianter Tensor 
Dg,5 = (g up Va O^ + SuaVp%") de (7) 


Das Symbol fy, bedeutet hier die kovariante Ableitung mit dem metrischen Tensor 
geg; über sich wiederholende griechische Indizes soll von 1 bis 4 summiert werden. 
Wenn das geometrische Obiekt & eine Funktion M geometrischer Obiekte y,, ist 


D = D (v, )) 
so drückt sich sein Liesches Differential durch die Formel 
DE = 21 Dy] P | (9) 
aus, wo i 
D 
br = -y A. Dyra | (10) 


als ein partielles Liesches Differential des Obiektes ® nach den Obiekten y,, benönnt 
wird. 3 | | 
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3. Das Variationsprinzip 


Aus der allgemeinen Relativitätstheorie ist bekannt, dass ein N ichtinertialbezugs- 


D 
system durch vier Funktionen £ (x”) charakterisiert werden kann. Aus ihnen bauen 
wir den symmetrischen metrischen Tensor 


BE JÉ 
Bab = = Ira Dat (11) 
auf. 
Im Inertialbezugssystem hat er die Gestalt 
Zap = (Eap)o: (12) 


wo (gap)o = 1 ist für a = 6 = 1, 2, 3, 4 und 0 für a + f. 

Betrachten wir jetzt eine Flüssigkeit, deren Zustand durch N Zustandsfunktionen 
die die Komponenten M geometrischer Obiekte y,, sind, beschrieben ist, Nehmen 
wir an, dass nur vier von diesen Funktionen unabhängig sind und dass die übrigen 
sich durch diese vier mittels integrabler oder nichtintegrabler Relationen ausdrücken 
lassen. Unsere Flüssigkeit werden wir ausserdem durch eine Wirkungsfunktion 


S= f DV— g (da), g= det (g,,), (13) 
charakterisieren, wo ® eine inwariante Funktion der AE gag und 
der Zustandsfunktionen 74 ist | 

D = Ó (bay; Vra): (14) 
Es bezeichne ó [gel (G V— a): die Variation von G / g, die der Variation der 


Inertialfeldpotentiale ó gag entspricht. Als Impuls-Energie- Tensor der betrachteten 
Fliissigkeit bezeichnen wir den durch die Formel 


è [gag] (© Y— 8) = TA Y— 588, (15) 


bestimmten ‘Tensor Te. Bei Beachtung, dass ó [g,] V— g=4V—ger ÒBap» 
erhalten wir aus (14) und (15) 


PWE DW a (16) 
9g 2 
Die Gleichungen , 
Va TB — 0 (17) 


heissen Bewegungsgleichungen. | pini 
Jetzt beweisen wir den folgenden Satz: Aus der Forderung, dass die Wirkungs- 


funktion (13) bei beliebiger Raumtransformation o gegenüber dem Lieschen 
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Partialdifferential D{y,4] S stationär ist und aus der Annahme, dass am Rande des 
Integrationsbereiches die Raumdeformation sowie ihre Ableitungen entsprechender 
Ordnung verschwinden, folgen die Bewegungsgleichungen (17). 

Tatsächlich folgt mit (9) aus (13) 


[Pieas: md 9 Y— g (dx) = f Dag] 9 V— g (da) + [Divi] 9 Y— (dx). (18) 
Aus (5) geht 
f Dis vral 9 Y — (di) = 0, (19) 
hervor und damit bekommen wir aus (18) 
DI V— a (dx) = — f Dia] P V— g (dy). (20) 
Bei Beachtung von (8) erhält man 
D [gas] 5 = i Dgap (dx). (21) 
was wir auch in der Form 
D [via] S=—2 f T” g Vaw" de (dx) (22) 


schreiben kénnen. Die partielle Integration führt zu 


Dy] S=—2 f ge p, T” Y— g (dx) de (23) 
Bei Beachtung von 
D[vral 5 = 0 (24) 
erhältman 
f aset, T” V— g de (dx) = 0. (25) 


Die Willkür des Vektors wf liefert die Bewegungsgleichungen (17). 


4. Bewegungsgleichungen für ideale und zähe Flüssigkeiten 


A. Ideale Flüssigkeit 


Eine ideale (reibungslose) Flüssigkeit kann durch folgende Funktionen beschrieben 
werden: a) kinematische Funktionen u“, die einen kontravarianten Vektor — den 
Geschwindigkeitsvektor des Medium — bilden; b) thermodynamische Funktionen 
E,W,S,F,T,p und n, die Skalare sind und entsprechend innere Energiedichte, 
Wärmefunktionsdichte, Entropiedichte, thermodynamisches Potential, Temperatur, 
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Driick und Dichte des Medium bedeuten (n ist die Dichte im Bezugssystem, in welchem 
das Volumenelement des Medium sich in Ruhe befindet — s. z. B. Landau und Lifszic 
1954, § 125). 


Die kinematischen Funktionen u? erfüllen die Gleichung 


Sail (26) 


Daraus folgt, dass nur drei Vektorkomponenten u? unabhängig sind. Unter den thermo- 
dynamischen Funktionen sind nur zwei voneinander unabhängig; die übrigen drücken 
sich durch diese zwei mittels der Zustandsgleichung und entsprechender thermody- 
namicher Identitäten aus. Daher haben wir fünf unabhängige Funktionen, zu deren 
Bestimmung wir ebensoviele unabhängige Gleichungen brauchen. Die fünfte setzen 
wir in der Form einer Kontinuitätsgleichung 


Pa (nu) = 0 


an (Falls diese Gleichung nicht erfüllt ist, so ist das thermodynamische Potential F gleich 
Null und dementsprechend haben wir dann nur vier unabhängige Funktionen zu 
bestimmen — s. z. B. Khalatnikov (1954)). 

Wir wollen jetzt zeigen, dass die bekannten Bewegungsgleichungen der idealen 
Flüssigkeit aus dem von uns formulierten Variationsprinzip folgen, wenn wir an- 


setzen 
S= f [W (Ve, u — 1) — p] Y— g (da) : (28) 
oder 
d — W (Ves u u. — 1) —p, E 
wo W= E +p. 


Fiir das Liesche Differential der Wirkungsfunktion (28) bekommen wir den Aus- 
druck 


d ID IDIN IT 
DIW, p, ws- | [sg DW + 5 Dp + pie Du IEC (30) 


Die im Integranden vorkommenden Ableitungen haben die Gestalt 


96 S Ne renn 31 
F En gemit ET (31) 


Durch Einsetzen von (31) in (30) erhalten wir bei Beachtung von (4), (5) und (6) 
DV, p, u] 5 = -f gag ul u" p pe V— g de (dx) + 


+f Ye per pape gag uf y u] wo V— g de (dx). (32) 
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2 TTT 


Wenn wir im ersten Glied auf der rechten Seite eine partielle Integration ausführen 
und (26) berücksichtigen, so erhalten vir 


I j: 284,0" Va T? V— gde (dx), | (33) 
wo 
Te. Pu ul — pg. (34) 


Das Verschwinden der Variation (33) und die Willkür der Deformation w“ führen 
zu den Bewegungsgleichungen 


. (Vu. pg. 0. (35) 


Eine einfache Rechnung zeigt, dass der Tensor T" mit dem Impuls-Energie-Tensor 
identifiziert werden muss. Die Gleichungen (26), (27), (35), die Zustandsgleichung 
und die thermodynamischen Identitäten bilden ein vollständiges Gleichungssystem 
fiir eine ideale Fliissigkeit. 


B. Zähe Fliissigkeit 


Für eine zähe Flüssigkeit nehmen wir die Wirkungsfunktion (13) in folgender 
Gestalt an 


S= f Rz uau? —1) — + 60 118 | V=8 (da), (36) 


wo der kontravariante Tensor II"? Spannungstensor heisst. In diesem Fall hat die 
Funktion ® die Form 


T i 
P= W (Vs us yes 1) 5 Sab Iet. | (37) 


Fiir den Impuls-Energie-Tensor bekommen wir nach (16) 


1 a 
Ted — en W Un Ub, ge naeh TI). (38) 


Die Bewegungsgleichungen lauten daher 


il 
A, (Wut ub — > Eu II” gaB — Nes) = 0. (39) 
Aus dem oben bewiesenen Satz geht hervor, dass : 
DIW, us, Defi S — 2 f &,, 0" y, TP V— de (dx). (40) 


Die Gleichungen (26) und (40) bilden kein vollständises System. Es miissen nach 
wectere Bedingungen hinzugefügt werden, die die Funktionen 7778 durch u? und die 
thermodynamischen Funktionen auszudriicken gestatten. Als zusätzliche Bedingungen 


PE tons 
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kann man z. B. analoge zu denen fiir den nichtrelativistischen Fall annehmen, die 
zu den Navier — Stokesschen Gleichungen fiihren (s. z. B. Kuzniecow, 1952, kop. 2.). 


Zum Schluss bemerken wir, dass das von uns angegebene Variationsprinzip 
seine Giiltigkeit auch im Gravitationsfeld bewahrt. 


Herrn Prof. R. S. Ingarden danke ich herzlich für wertvolle Ratschläge bei dieser 
Arbeit, sowie für das fòrdernde Interesse, das er ihr entgegenbrachte. 
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THEORY OF ORDER IN Fe,Al ALLOY 


KAZIMIERZ F. WOJCIECHOWSKI 


Department of Experimental Physics of the Wrocław University 
(Received April 2, 1958) 


A statistical theory of order — disorder transitions in the F ez Al alloy is presented 
which enables us to take into account the interactions of atoms from farther neighbourhood 
and thus to determine more accurately the transition temperature from order to disorder. 
It is shown that the interaction energy of third neighbourhood is approximately 1/2 of 
the interaction energy of the second one. 


1. Introduction 


Although the theory of binary alloys has been widely developed (Muto et al. 
1955), the problem of the precise determination of the transition temperature from 
order to disorder remains still open. Besides, there were till now no theoretical in- 
vestigations concerning especially the Fe,Al alloy which shows many interesting 
features connected with orderly arrangement and becomes ferromagnetic in a certain 
temperature interval. The main aim of this paper will be therefore to present a statistical 
theory of order — disorder transitions in this alloy. This theory will be a generalization 
of the theory of Bragg and Williams (1934) and will enable us to take into account the 
interactions of atoms from farther neighbourhoods than the first. 


2. Fundamental Equations and Long-range Order Parameters 


If we denote by a,, the number of atoms of the i-th kind in the lattice points of the 
u-th sublattice, and m N the total number of atoms in the alloy, we receive for the 


FeAl alloy the following equations 


3 
du + 412 + Ms = a N. 


N 
422 e (1) 


Aj + lap = lig + laz» 
N 
dy = const = roy n 0. (2) 
(15) 


y 
á 
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In these, i = 1 corresponds to Fe atoms, i = 2 to Al atoms and 4 = 1, 2, 3 denote 
the three sublattices into which we can decompose the superlattice of the Fe,Al 
alloy (Fig. 1). Relations (2) hold because the transition from order to disorder concerns 
only atoms placed-in the lattice points of the sublattices 2 and 3. 


Fig. 1. Superlattice of the FeAl alloy and its decomposition in thre esublattices. o = 1 = a Va, 07 
=2=a,0=3=a y2 


On the basis of the general definition of the long-range order parameter given by 
the author (1958), we shall define this parameter in the present case as follows 


TL _ : " (3) 


For s — 1, there is complete order in the alloy, for s — 0, complete disorder. From 
equations (1) and (3) we shall express all a;, as functions of s 

Ay = dp = Ea +s), 
av = am = (1 — 9. (4) 
Configurational Energy 


In the paper referred to above the author has given a general expression for the 
configurational energy. In the case considered here this formula can be written in the 


form 
2 3 i 
2 2 * yo Y wo aru pj Vig (Qs ` | (5) 


e ki py 
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where o — 1, 2, 3 denote the distances between first neighbours, second neighbours 
etc. respectively, 2,,(Q) is the number of nodal points of the v-th sublattice for which 
the distance from a given point of the w-th sublattice is 0; qy; (0) is the number 
of pairs consisting of kj atoms, which mutul distance is 0; Py, is the probability 
that a lattice point of the H- th sublattice is occupied by the atom of the k-th kind, 

Pku = 75 


u 


(6) 
where M, is the number of nodal points of the u-th sublattice, V,; (p) is the interaction 
energy of atom k with atom j, their distance being o. 

Formula (5) enables us to calculate E, to different orders of approximation. In 
our case we shall calculate E taking into account the interactions of atoms from the ' 
third neighbourhood (o = 3). In order to do this, let us first calculate E (1) i. e. the 
configuration energy in which only interactions of the nearest neighbours are taken 
into consideration (o = 1). From the Fig. 1 it is clear that 

212 (l) = 213 (1) = 4, 
Z (1) = za (1) = 8, 
zu (1) = 255 (1) = 333 (1) = Zas (1) = za (1) = 0 (7) 
Making use of formula (5), we find 
E (1) = const. (8) 
Since for 9 = 2 (Fig. 1) 
Z(2) = zu (2) = sás (2) = 232 (2) = 6, 


212 (2) = 232 (2) = 212 (2) = 221 (2) = 222 (2) = 233 (2) = 0 (9) 
we receive 
3 
1 1 2 AE 
qu (2) = 5 2% (2) aku Pir = > [Z11(2)arı Dj1 + 223 (2) % Piz + 282 (2) akz Djs 
= 3 (ay Pa + d Pija + Urs Djs) (10) 


Thus 
3N 
Ou (2) = FI (5 — $2), 


ou © = 22 A— 5%), 
du (2) = A (1 + só). 


‘Substituting these values in (5), we obtain 


EO = rate + Val HIM EL AY 


y 
7 
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ÄR Tele 
where V (2) is expressed in terms of the following relation 


V (0) = 27/1 (0) — Pile) — V220), (12) 
foro. 
For o = 3, since 


Z (3) = zy G) = ze (3) = za (3) = 12 
223 (3) = 232 (3) = 213 (3) = zz (3) = 212 (3) = za (3) = 9, 


we obtain analogously 


3N 

an @) = E 6 +59, 
3 

ou (8) — a + 


mB — a — 9), 
and 
EB = — FN (5 Vn 8) + Vaa B) +2 V G- Fiat, ` W) 


where V (3) is given by (12) for o = 3. 
Thus, taking into account the interaction of the third neighbourhood we get finally 


EGO = E(1) + £2) + £3) = — SN (V(2) — 2V(3)) s? + const. (14) 


4. The State of Thermodynamical Equilibrium and the Determi- 
nation of T,. 


The state of thermodynamical equilibrium is defined by the condition of minimum 


configurational free energy F (s) = E (s) — kT In W (s): 


OF (s) _ 
5 (15) 


where W (s) is the thermodynamical probability. Using the value W (s) calculated 


according to the Bragg- Williams approximation for the AB-alloy and relation (14), 
we receive from (15) 


=s 2 kT 


This equation corresponds to the equation received by Bragg and Williams (1934) 
for the AB-alloy, but thanks to the taking into account of the third neighbourhood 
interactions, it enables us to determine the transition temperature from order to disor- 


do adi SE TASK ae 
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der T, more accurately. The temperature T, can be obtained from (16) in the follow- 
ing manner: 


lim (2 niti)- E 
e 1-5 IKT, 


s- 0 


hence 


31. — 2V(3)] 


Leo. Ak 


(17) 


5. Comparison of the Theoretical Results with Experiment 


The Fe, Al alloy was the subject of many experimental investigations. The transition 
temperature from order to disorder for the Fe,Al alloy found in the most recent of 
them (Rode 1957) is 530°C. Assuming that the taking into account of the third neigh- 
bourhood interaction enables us to determine correctly T, by means of (17), we can 
assume that 


3 530°C k = V(2) — 2V(3) 


and hence calculate the difference of the reduced energies of interaction between the 
second and third neighbourhoods. Taking for k = 1.380 x 10-16 g cm? sec”? deg`!, 
we find that 

V (2) — 2V (3) = 0.0058 eV 


hence 


V(3) = = V (2) — 0.0029 eV 


` Thus approximately 
1 
ve) K s FO), 


i. e. the interaction energy of the third neighboruhood can perform a considerable 
influence on the ordering processes and should not be neglected.` 
Although the theory presented here enables us to determine 7, more precisely, 


the agreement with the experiment is only a qualitative one. Fig. 2 shows for example 


T 


the experimental curve and the theoretical dependence s = f Tl 


It is worth while to remark that on the basis of the experimental results of Rode 
(1957) one can determine the constant A for the Fe; Al alloy which appears in formula 


(25) of Smirnov's paper (1947), Now, since near of T, 
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Fig. en ot e Sch =f (27. Continuous ling — experimental en of Rode (1957) E 
broken line — theoretical curve. l b 

the specific resistance o (T) will be expressed on the basis of fortu Co) of Smirnov’s 

paper (1947) as follows | té AO Dee 

$ «n- (5-7 : 

In the paper of Rode o (T,) = 125 u Q em, thus Ape, 4; = 500 u Q cm. 

At last we want to call attention to the fact that all the quantities obtained Ee 

do not depend upon the interactions of the first neighbourhood. This fact is in agree- . 

ment with experiment, since it has been established that only MEC in the 

latice pani of the sublattices 2 and 3 take part in the process of or g 


Spire . Wille e E. au 
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UBER DEN ZUSAMMENHANG DER 
BEWEGUNGSGLEICHUNGEN MIT DEN FELDGLEICHUNGEN 


VON CZESLAW JANKIEWICZ 


Fakultät fär Theoretische Physik, Universität, Wroclaw: 
Mathematisches Institut der Polnischen Akademie der Wissenschaften, Wroclaw. 


(Eingegangen am. 26. Mai 1958) 


Bei Verwendung des Lieschen Differentiales wird gezeigt, dass in jeder klassischen 
auf Grund des allgemeinen Relativitätsprinzipes formulierten Feldtheorie die Bewegungs- 
gleichungen der kontinuierlich verteilten Feldquellen aus den Feldgleichungen folgen und 
dass die Bewegungsgleichungen mit den Erhaltungssätzen des Feldes identisch sind. Als 
Beispiel dient die Gravitationstheorie und die: Elektrodynamik. ) 


1. ‚Einleitung 


In der klassischen Feldtheorie unterscheidet man im allgemeinen zwei Funktions- 
arten, die den Zustand des physikalischen Systems beschreiben. Funktionen, die die 
Bewegung der Quellen des betrachteten Feldes charakterisieren — wir werden sie 
Bewegungsfunktionen nennen — und Feldfunktionen, die das Feld der betrachteten 
Quellen charakterisieren. Dabei gibt es zwei Betrachtungsarten der Feldquellen. 
Die erste beruht darauf, dass man die Quellen als Unstetigkeiten und Singularitäten 
. des Feldes ansieht, die zweite stützt sich auf die Annahme der stetigen und regulären 
Verteilung der Feldquellen. Die erste Betrachtungsart kann man als speziellen Fali der 
zweiten ansehen, falls man Quellenverteilungen vom 6- Typus zulässt. Daher 
werden wir im Weiteren die F eldquellen ausschliesslich als stetig ‚verteilt 


i 


beschreiben. | 

Die Ableitung der Gleichungen fär die Bewegungs- und Feldfunktionen aus dem 
Experiment ist eine der wichtigsten Aufgaben der Physik. Die Zahl der Gleichungen 
muss eine solche sein, dass man aus ihnen bei gegebenen zusätzlichen Bedingungen 
(Anfangs- und Randbedingungen) eine eindeutige Lösung zu erhält. Da im allge- 
meinen die Bewegung der Feldquellen von den Veränderungen des F eldes, und die 
Feldveränderungen von der Quellenbewegung abhängen, bilden die Gleichungen 
für die Bewegungs- und Feldfunktionen ein Gleichungssystem mit unseparierbaren 
unbekannten Funktionen. Daher ist die Aufteilung der Gleichungen in Feld- und 
Bewegungsgleichungen nicht eindeutig. Über einen Zusammenhang zwischen den 


(21) 
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Bewegungs- und Feldgleichungen kann man erst dann sprechen, wenn man sie vorher 
definiert hat. Speziell hat z. B. die Feststellung, dass fiir ein physikalisches System die 
Bewegungsgleichungen aus den Feldgleichungen folgen und fiir ein anderes System 
das nicht der Fall ist, nur dann einen Sinn, falls in beiden Fallen die Aufteilung der 
Gieichungen nach denselben Kriterien durchgeführt wurde. 

Einstein und Grommer (1927) haben zum erstenmal gezeigt, dass in der allgemeinen 
Relativitätstheorie die Bewegungsgleichungen eines Probeteilchens, d. h. die 
Gleichungen einer geodätischen Linie aus den Gravitationsfeldgleichungen folgen. 
Weiter haben Einstein, Infeld und Hoffmann (1938) die Bewegungsgleichungen fiir 
ein System materieller Körper aus den Feldgleichungen hergeleitet, und zwar in 
einer Naherung, die über die Newtonsche Theorie hinaus geht. Unabhangig von ihnen 
hat auch Fock (1939) auf Grund einer anderen Methode die Newtonschen Bewegungs- 
gleichungen erhalten. Petrova (1949) und Papapetrou (1951) haben mit dieser Methode 
höhere Näherungen berücksichtigt. Eine andere Form der Bewegungsgleichungen 
haben nachher Einstein und Infeld (1949) abgeleitet. Ihre Rechungen waren sehr 
schwierig und verwickelt und erst Infeld (1954) hat sie durch die Verwendung der 
6-Funktion sehr vereinfacht. 

Die Identitäten von Noether (1918) spielen eine grundsätzliche Rolle bei der 
Ableitung der Bewegungsgleichungen aus den Gravitationsfeldgleichungen. Bergmann 
(1949) und Goldberg (1953) haben gezeigt, dass in der Gravitationstheorie diese Identi- 
täten vom metrischen Tensor unabhängig sind und dass in einer jeden kovarianten 
Feldtheorie die Bewegung der Feldsingularitäten aus den Feldgleichungen bestimmt 
werden kann. 

Infeld und Wallace (1940) haben gezeigt, dass in der klassischen Elektrodynamik 
die Bewegungsgleichungen aus den Maxwellschen Gleichungen folgen; jedoch ist 
die Ableitung nicht eindeutig. Dieser Mangel an Eindeutigkeit beruht auf der Willkiir- 
lichkeit von vier Funktionen, deren Auswahl erst die Herleitung der eee 
gleichungen ermöglicht. 

In der vorliegenden Arbeit wird das allgemeine Relativitätsprinzip formuliert ohne 
zuhilfenahme der für die Gravitationstheorie spezifischen Begriffe. Im Zusammenhang 
damit werden an Hand eines Variationsprinzips kovariante Bewegungsgleichungen 
für die kontinuierlichen Quellen definiert. Es wird gezeigt, dass in jeder klassischen, 
auf Grund des allgemeinen Relativtatsprinzipes formulierten Feldtheorie diese Be- 
wegungsgleichungen aus den Feldgleichungen folgen. Es wird bewiesen, dass die 
Erhaltungssätze, die aus der Kovarianz der Feldgleichungen folgen, identisch mit den 
Bewegungsgleichungen sind. Diese Tatsachen werden am Beispiel der Gravitations- 
theorie und der Elektrodynamik erlaütert. Speziell wird gezeigt, dass das allgemeine 
Relativitätsprinzip die Gestalt beliebiger vier Funktionen in der auf die Elektrodynamik 
angewendeten Infeld-Wallace-Methode bestimmt. 
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2. Das Liesche Differential 


Als Unterlage fär unsere geometrischen Betrachtungen soll ein Puktraum X, dienen 
(Golab 1956). Eine der grundsätzlichsten Fi i e 

| a r grundsätzlichsten Eigenschaften eines solchen Raumes ist 
die, dass man jeden seiner Punkte einen Punkt im n-dimensionalen Raum der 
Zahlen xt, x?,..., x” zuordnen kann. (Die durch griechische Buchstaben bezeichnten 
Indizes nehmen die Werte 1, 2,..., n an). Diese Zahlen nennen wir die Punktkoordi- 
naten des Raumes X,. Die Zuordnung neuer Zahlen à den Punkten a* des 
Raumes X, gemäss den Formeln 


, : O 
* — x (x*), det = ) + 0 (1) 
heisst eine Koordinatentransformation. 

Mit dieser Transformation ist der Begriff des geometrischen Obiektes verbunden 
(Wundheiler 1934). Ein System von N Funktionen %%) (die durch grosse lateinische 
Buchstaben bezeichneten Indizes durchlaufen die Zahlen 1, 2,..., N) bildet die Kom- 
ponenten eines geometrischen Obiektes, wenn folgende Bedingungen erfiillt sind: 

1. Nach der Koordinatentransformation (1) sind die neuen Komponenten des 
geometrischen Obiektes Funktionen der alten Komponenten dieses Obiektes, der 
alten Punktkoordinaten, der neuen Punktkoordinaten sowie ihrer Partialableitungen 
nach den alten. Das Transformationsgesetz fiir das Obiekt schreiben wir kurz in der 
Form 


Dig (aes |) == OA Wa IA); 
2. Die Funktionen G, erfiillen die Gruppenbedingungen 


CAG (vs net) ott | == CP ap Ga) 
Ga [Gar (Pa, x*, x), *, a*] = Pa 
G 4 (Wa: *, X7) = Ya: (3) 
Die geometrischen Obiekte heissen linear, wenn die Transformationsgesetze die 
Form haben 


Pa 22, GyB (, x) vg + Ga (x^, SEH x (4) 
B 


wie z. B. bei Vektoren, Tensoren, Tensordichten und Christoffelschen Symbolen. 

Neben den Koordinatentransformationen (1) kann man auch einparametrige 
Punkttransformationen des Raumes X, oder die ihnen entsprechenden infinitesimalen 
Transformationen 


* — x FW" de (5) 


betrachten, wo w“ (x”) ein kontravarianter Vektor und de eine infinitesimale Grösse . 
ist. Mit diesen Transformationen steht der Begriff des Lieschen Differentiales in 
Verbindung. der zuerst von Slebodzinski (1931) eingeführt wurde. Unter dem Lieschen 


24 i Czestaw Jankiewicz 


Differential eines geometrischen Obiektes verstehen wir den in de linearen Ausdruck 


D y4 = va (0) — va (), (6) 
wo 
"Wall = Ga [va (x“), x, x°]. (7) 
Daraus erhalten wir (s.z. B. Yano 1955) 


Dy e= (ot 97 = 0,070 srt Af 9, 9,0 — ...) de, (8) 


wo G 4» Gao» C me , +.. Patrialableitungen entsprechenden Grades der F a G4 
nach "xt, 89,'x?, 9 SEN Yon im Punkte de = 0 bedeuten; hierbei ist 9, = 9/9x%, 9, 9, = 
= 931927 D . (Über Ban wiederholende aka Indizes soll von 1 bis n 
summiert -— 

Bei linearen Obiekten hängen die Funktionen G4,, Ga", Ga"... linear von den 
Komponenten y4 ab. So drückt sich z. B. das Liesche Differential einer Tensordichte 
vom Gewicht +% durch die Formel aus 


DB," (we 9, B," —B "93 o*—pB 9, oi + 
+ By” o, e + B dw + k B, 9,09) de. (9) 
Eine Verallgemeinerung des Begriffes eines geometrischen Obiektes ist die geo- 
metriche Komitante (Gołąb 1938). Betrachten wir eine Funktion G, die von einem 


geometrischen Obiekt y, und seinen Ableitungen nach x* abhängt. Sie habe im jeden 
Koordinatensystem die Form 


D =P (pa, 9, Va» Ia Ig Vas). (10) 
Wenn bei einer Transformation des Obiektes y, die Funktion & sich ebenso wie die 
komponente eines geometrischen Obiektes transformiert, so heisst sie eine geometri- 
sche Komitante. (Wir beschränken uns auf einkomponentige Komitante, die für 
unsere physikalische Betrachtungen ausreichen). Man kann zeigen (Jankiewicz 1958), 
dass das Liesche Differential der Komitante (10) sich durch die Formel 


fico Di JD : 
Do = SI e rp di 
Bundes lasst. 


; Wenn man den Ausdruck (11) beiderseits úber ein Bereich des Raumes ON 
integriert und dabei annimmt, dass sowohl der Vektor w” als auch seine ible tigen 


entsprechenden Grades nach x* am Rande des Integrationsbereiches verschwinden, 
so erhält man 


il DÖ (dx) = f ba L4 DY , (da), | (12) 


wo 


S 96 96 96 
ai Cáisc ah ei WO u an ae ape n 
ee ee TC E eg dà 


W 
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als funktionale oder Lagrangesche Ableitung bezeichnet wird. Wenn wir (8) in Be- 


tracht ziehen und auf der rechten Seite von (12) eine partielle Integration durchfiihren, 
so erhalten wir 


f 0949 = f (2,7 — Eier de (da), (14) 
wo 
D SS = [22 e m 92 (La Ca SE d, de (La ep) car e] (15) 
A 
DE 2 LA (CA. — 9, YA). co 


Jetzt hänge die geometrische Komitante ® von M geometrischen Obiekten und 
ihren Ableitungen ab: 


d — O (l- 0 A- 9 Og ps). (17) 


(Von nun an sollen der Index I Obiekte, dagegen der Index 4 Komponenten 
numerieren, wobei sie entsprechend verschiedene Werte durchläufen). Das Liesche 
Differential der Komitante (17) drückt sich durch die Formel 


DP 2 D[vialP (18) 


aus, wo D[y,,1® ein Liesches Partialdifferential ist, das für jeden Wert von I durch 
(11) gegeben ist (Jankiewiez 1958). Aus (12) und (18) erhalten wir 


f DG (dx) = f DD Lra Dyra (dè. (19) 
TUA 
Führen wir auf der rechten Seite eine partielle Integration aus, so bekommen wir 
f DG (dx) = f X (9, Tr" — Fio! de (da). (20) 
T 


Die Ausdrücke L,,, Tp“ und Fi sind für jeden Wert von I entsprechend durch die 


Formeln (13), (15) und (16) gegeben. 


3. Bewegungsgleichungen und Feldgleichungen 


Geometrisch ausgedrückt, können wir das af!gemeine Relativitätsprinzip (Ein- 


stein 1916) folgendermassen formulieren: 


a. Das physikalische Raum-Zeit-Kontinuum ist ein vierdimensioneller Raum Xo 
wobei x}, x2, x8 die räumlichen, x* = it die zeitliche Koordinate bedeutet (i = y: 

b. Die Raum-Zeit-Koordinaten zweier (im allgemeinen nicht inertialer) Bezugs- 
systeme sind mittels der Formeln (1) verkniipft. 

c. Die physikalischen Gròssen sind geometrische Obiekte. | ! 

d. In allen Bezugssystemen haben die physikalischen Gleichungen, die den 


Zusammenhang zwischen physikalischen Grössen ausdrücken, dieselbe Gestalt. 
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Falls wir in den obigen Annahmen die Transformationen (1) durch Lorentztrans- 
formationen ersetzen und uns auf Inertialbezugssysteme beschranken, so erhalten 
wir das spezielle Relativitätsprinzip. 

Im Zusammenhang mit dem allgemeinen Relativitätsprinzip nehmen wir an, dass 
die Bewegungsfunktionen Komponenten eines geometrischen Obiektes y sind 
und die Feldfunktionen — Komponenten M geometrischer Obiekte “ra (x). 

Das Wirkungsfunktional des betrachteten physikalischen Systems postulieren 
wir in der Form 


8 f D (dx), (21) 


D =D (wg, Vra» Va VIA» da VIA») (22) 


eine geometrische Komitante ist, die von Bewegungs- und Feldfunktionen sowie den 
Ableitungen der letzteren nach x* abhängt. Wir betonen, dass in (22) die Ableitungen 
der Bewegungsfunktionen nicht auftreten. 

Wir definieren die Bewegungsgleichungen als Gleichungen, die aus der Forderung 
der Stationarität des Wirkungsfunktionals D[y,1S = 0 folgen mit dem Ansatz, dass 
am Rande des Intergrationsbereiches der Vektor w“ und seine Ableitungen nach & 
verschwinden. Bei Beachtung von (14) bekommen wir fir (21) 


Diy,] S = P (8, T — Vi o de (da). (23) 
Infolge der Willkür von w” erhalten wir die Bewegungsgleichungen 
et, (24) 


Das Obiekt T., werden wir Impuls-Energieobiekt und F, — Kraftobiekt der Feld- 
quellen nennen. Die Gleichungen (24) kann man als (verallgemeinertes) Erhaltungs- 
satz fur Impuls und Energie der Feldquellen interpretieren. Im Einklang mit diesen 
Definitionen sind also die Bewegungsgleichungen identisch mit den Impuls-Energie- 
Erhaltungssätzen der Feldquellen. 

l Dementsprechend führt die Forderung der Stationarität des Wirkungsfunktionals 
(21) gegenüber der Variation der Feldfunktionen öy,, und die Annahme, dass am 
Rande des Integrationsgebietes diese Variationen sowie ihre Ableitungen verschwinden, 
direkt zu den Feldgleichungen. Aus (21) folgt 


o[vr4] S = f 2 2 La övra (da). (25) 
Die Feldgleichungen haben also die Gitali | 


Lia =0 (26) 


Die Ansätze des allgemeinen Relativitätsprinzipes werden erfüllt, wenn ® eine 
Skalardichte mit dem Gewicht +1 und y, und yr, linear sind. Diese Bedingungen 
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sind ausreichend, da dann die linken Seiten der Gleichungen (24) und (26) homo- 
gene lineare Obiekte sind. 
Bei Beachtung von (12) und (18) folgt für (21) 


f Diva vral P (da) = (IE. T^ — Harde + DD Lia Dv] (d). (27) 
I. A 


Aus (9) geht hervor, dass das Integral eines Lieschen Differential einer Skalardichte 
vom Gewicht + 1 identisch verschwindet 


f Diva Yra] (dx) =0 (28) 
und das führt (27) in 
/ [(2,, 7. — E)e'de +3 Y Lj Dyral (dx) =0 (29) 
1 A 


über. Wenn die Feldgleichungen (26) erfüllt sind, dann liefert (29) die Bewegungs- 

gleichungen (24). Auf diese Weise haben wir einen allgemeinen Beweis dafür erhalten, 

dass die Bewegungsgleichungen aus den Feldgleichungen folgen. så 
Mit (20) nimmt die Identität (29) die Gestalt 


J (e, T^ —F) +2 (9, 2h — Fó] a" de (dx) zo (30) 


an, und daraus ergeben sich die Identitäten von Noether (1918) in einer zur physika- 


lischen Interpretation besonders geeigneten Form 


(% 7. — F) + D (% Té — Fr) = 0 (31) 


(Diese Identitäten stellen eine Verallgemeinerung der strengen Erhaltungssätze (Berg- 
mann 1949, Goldberg 1953) auf Felder und Feldquellen zugleich dar, die mithin nicht 
unbedingt durch differentielle geometrische Obiekte der klasse 1 (Yano 1955) be- 


schreiben werden miissen). | 
Die Obiekte T,“ werden wir Impuls-Energieobiekte und Ii, — Kraftobiekte. des 


Feldes nennen. Die Gleichungen 
2 (0% Tj," — Fr) = 0 (32) 
kann man als (verallgemeinerte) Erhaltungssätze für Feldimpuls und Feldenergie 
interpretieren. Aus (31) folgt ein allgemeiner Beweis für die Identität der Bewegungs- 
gleichungen (24) mit den Erhaltungssätzen für Feldimpuls und Feldenergie (32). 
Diese Identitäten gelten auch dann, wenn die Feldgleichungen (26) nicht erfüllt sind. 
Es sei ein physikalisches System durch n Bewegungs funktionen y, und ny + na + 
+... + ny Feldfunktionen A charakterisiert. Für diese Funktionen haben wir 
n Bewegungsgleichungen (24) und n, ++... TM Feldgleichungen (26). 
Wie wir gezeigt haben sind die Bewegungsgleichungen nicht unabhängig sondern 
folgen aus den Feldgleichungen. Daher haben wir bloss n4 + n5 4- ... + ny, unabhängige 
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Gleichungen für ni m+n E. E ny Funktionen, d.h. die Zahl der unbekannten 
Funktionen ist um n gròsser als die Zahl der unabhängigen Gleichungen. Um die 
Lösung des Problems eindeutig zu machen, muss man zu den Feld — und Bewegungsg- 
leichungen zusätzliche Bedingungen hinzufügen. Vom physikalischen Standpunkt aus 
bedeutet die Auslesung dieser Bedingungen die Wahl eines Bezugssystems. Es werden 
zwei Methoden zur Wahl des Bezugssystems verwendet. Die erste von ihnen beruht 
darauf, dass man n zusätzliche Gleichungen für die Feldfunktionen (bzw. direkt n 
Feldfunktionen) und Anfangs — und Randbedingungen, die die eindeutige Lösung 
aller Gleichungen garantieren, angibt. Dieses Verfahren wurde von Fock (1939) in 
der Gravitationstheorie angewendet und — wie wir zeigen werden — von Infeld und 
Wallace (1940) in der Elektrodynamik. Die zweite Methode fusst auf der Anwendung 
eines entsprechenden Näherungsverfahrens, das die Eindeutigkeit der Lösung der 
Feld- und Bewegungsgleichungen sichert. Sie wurde von Einstein, Infeld und Hoff- 
mann (1938).in der Gravitationstheorie angewendet. 


4. Die Gleichungen des Gravitationsfeldes 


In der allgemeinen Relativitätstheorie bildet das Trägheits- und das Gravitations- 
feld zusammen das verallgemeinerte Gravitationsfeld. Das Feld ist durch den symme- 
trischen metrischen Tensor ga, (x”) eines vierdimensionellen Riemannschen Raumes 
beschrieben. Das Bogenelement dieses Raumes hat die Form 


ds? = ga, da® dé. (33) 


Wir nehmen an, dass die Bewegung der Quellen des verallgemeinerten Gravitations- 
feldes durch eine kontravariante Vektor — dichte p* vom Gewicht + 1 charakterisiert 
ist. Die Grösse p = V ga, p* p? stellt die Quellendichte, und u“ = o^! p* — die Quellen- 
geschwindigkeit dar. Der Vektor u“ erfiillt also die Identitat 

gag U^ uf = 1. (34) 


Demnach stellen die Tensorkomponenten g,, die Feldfunktionen dar und die 
Komponenten der Vektordichte p* die Bewegungsfunktionen. 
Das Wirkungsfunktional nehmen wir in der Gestalt 


Se f E Sap p* PP — 1225 SR ZE | (da), 65 


an, wo R = g,, R“ die Invariante des Ricci — Tensors ist; c, m und y sind physi- 
kalische Konstanten (Lichtgeschwindigkeit, Masse und Gravitationskonstante). 

Das Liesche Partialdifferential nach der Bewegungsfunktion p* ist für (35) durch 
die Formel 


D[p*] 5 = f mg, u* D pP (dx) p (36) 
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gegeben. Bei Beachtung von (9) und (34) 


ben. bekommen wir nach einer partiellen Inte- 
gration in der letzten Formel 


D [p°] S = [me ga (ui p, u* N ul p, p^) w de (dx), (37) 


wo f7, die kovariante Ableitung mit dem metrischen Tensor 848 bezeichnet. Aus (37) 
folgen die Bewegungsgleichungen 


mo ga, (ut put rs Va p") = 0. (38) 
Wenn man (34) kovariant differenziert erhält man die Identität 
Sap U^ Vy u = 0. (39) 
Die Relationen (38) und (39) liefern die Bewegungsgleichungen in der Form 
m o gar 1% = 0, p,p*=0. 40) 
wo 
we =u" p, us (41) 


die Beschleunigung der Feldquellen ist. 
Die Variation von (35) nach den Feldfunktionen Zap hat die Gestalt 


ó [Zag] S = a us uf NT re lar ab IE (dx) (42) 
Saß -— 9 d 16 ay 8 2 8 Sag 


Diese Formel gibt uns die Feldgleichungen 


us uP, (43) 


Reb Rev = — 


Wegen (9) kònnen wir das Liesche Partialdifferential nach den Feldfunktionen 
gag von (35) folgendermassen aufschreiben 


Es 2 : 
D [ess] S= fr ET ut uf iam — 8 88 (ee 
5 Rest) | wrde (dx). | (44) 


Daraus ist ersichtlich, das die Impuls-Energietensordichte des Feldes durch den 
Ausdruck 


= V-g j am d 45 
„ ges (Re 2 R& (45) 
gegeben wird und ' 
1 9 gap 
a Tab (46) 
TREE 


das Kraftobiekt des Feldes darstellt. 
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Aus (32), (45) und (46) bekommt man die Erhaltungssätze fär Impuls und Energie 


des Feldes in ersichtlich kovarianter Form 
p, T, = 0. (47) 


Das zweite Glied im Integral (35) ist eine Skalardichte vom Gewicht + 1, die nur 
vom Tensor ga, und seinen Ableitungen abhängt. Das Liesche Differential dieses 
Gliedes sicht folgendermassen aus 


D [as] f K N E (dx) = | 2 (m -3 Reet) & ot de (da), — (48) 


woraus ohne weiteres die Identitäten von Bianchi abzulesen sind 


Va ( 2 = Res) = 0. (49) 


(Die aus dem ersten Glied von (35) folgenden Identitäten sind trivial.) 

Die Berechnung der kovarianten Ableitung der beiden Seiten der Feldgleichungen 
(45) führt bei Beachtung von (49) zu den Bewegungsgleichungen (40). Anderseits kann 
man die Bewegungsgleichungen auch aus (47) und (49) ableiten. Wir sehen also, 
dass aus den Gleichungen des verallgemeinerten Gravitationsfeldes die Bewegungs- 
gleichungen folgen und dass die Identität der Impuls-Energie-Erhaltungssätze des 
verallgemeinerten Gravitationsfeldes mit dem Bewegungsgleichungen auch dann noch 
gilt, wenn die Feldgleichungen nicht erfiillt sind. 


5. Die Gleichungen des elektromagnetischen Feldes 


Das elektromagnetische Feld kann mit Hilfe eines kovarianten Vektorpotentials 
beschrieben werden, jedoch sind wir nicht im Stande mit ihm allein das 
allgemeine Ralativititsprinzip zu erfiillen. Tatsächlich kann man aus einem- Vektor 
und seinen Ableitungen nach x* keine geometrische Komitante aufbauen, die eine 
Skalardichte vom Gewicht + 1 wäre. Vom physikalischen Standpunkt aus ist das 
klar, da doch ein im nichtinertialen Bezugssystem betrachtetes physikalisches System 
neben den Potentialen des elektromagnetischen Feldes noch durch die Potentiale 
des Tragheitsfeldes charakterisiert werden muss. Um die Bedingungen des allgemeinen 
Ralativitätsprinzipes zu befriedigen, werden wir als ein verallgemeinertes elektro- 
magnetisches Feld ein Feld ansehen, das durch Potentiale des elektromagnetischen 
und des Trägheitsfeldes beschrieben wird. Jedoch unterscheidet sich das verallge- 
meinerte elektromagnetische Feld grundsätzlich vom verallgemeinerten Gravitations- 
feld. Im letzteren können wir wegen der Identität der Trägheits — und der Gravita- 
tionsmasse die Potentiale des Gravitationsfeldes nicht von denen des Trägheitsfeldes 
unterscheiden. Im ersteren dagegen sind die Potentiale des elektromagnetischen - 
Feldes deutlich getrennt von denen des Trägheitsfeldes. Mathematisch sieht das so 
aus, dass die Potentiale des verallgemeinerten Gravitationsfeldes ein geometrisches 
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Obiekt bilden (eine physikalische Grösse), die des verallgemeinerten elektromagne- 
tischen Feldes dagegen — zwei Obiekte (zwei physikalische Grössen). 
In diesem Fall sind also die Feldfunktionen die Komponenten des kovarianten 


Vektors A, () und die vier skalaren Funktionen £ (x”). Die ersten beschreiben das 
elektromagnetische Feld und die zweiten das Trägheitsfeld. 


D 
Aus den &(x”) bilden wir einen symmetrischen Tensor zweiten Grades 


353 
Des dé 
han = d xe IxB o (50) 


i= 


Wir setzen an, dass die Potentiale des Tragheitsfeldes in einem Inertialbezugssystem 
folgende Bedingungen erfiillen 


Yap = TE (51) 
wo (Yeslo = l für a = f = 1, 2, 3, 4 und 0 für a Æ B ist. 


Die Bewegungsfunktionen sind die Komponenten der kontravarianten Vektor- 


dichte p* vom Gewicht +1. Die Grösse o = Vy, p« p? stellt die Quellendichte dar 
und u? = o^! p* die Quellengeschwindigkeit. Der Vektor u? unterliegt der Identität 


Vag U^ uf = 1. (52) 


Das Wirkungsfunktional postulieren wir in der Form 


1 , 
SE f [sar — $ Aart ib Hes lis Hy 69 
wo : 
y= det (Yap) HE = V. — J. A,» Haß = Et pi TTR 


und e die elektrische Ladung bedeutet. Das Symbol p, bezeichnet die kovariante 
Ableitung mit dem Tensor y,g- * 

Das Liesche Partialdifferential von (53) nach den Bewegungsfunktionen p* drückt 
sich folgendermassen aus 


piss — [| nais - 2 As | Dr? (aa) (54) 
Mit (9) erhalten wir nach partieller Integration 
e * 
D|p*] 5 = il [ne Ve (u^ Vn u* +u" Pap") — 0 — He o de (dx). (55) 
Aus (52), (55) und der Antisymmetrie des Tensors H,, folgen die Bwegungsgleichungen 


MOY av WE — 0 £ e (56) 
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Die Variation des Wirkungsfunktionals (53) nach den Feldfunktionen A,, € hat 
die Form 


i 1 E e 
Ô Aa, 95 V» He Var du Jaa, = 


1 
e E (moutw 4- V — y Uw) e| (dx), (57) 
wo 
pest Ss Vag H” Hb + i yi Hag Heb (58 ) 
An = 167% sé 


der Impuls-Energietensor des (nichtverallgemeinerten) elektromagnetischen Feldes 
ist. Bei Beachtung von (50) bekommen wir nach partieller Integration 


4 
i ie 
6 A, 8] S = Ne V—y y, H= 0 E 2 04,— 3: V» | mov w + 


i=1 


9x" 


Aus der Definition der kovarianten Ableitung geht hervor, dass 


+= um À | s] (da). 5 


a 
Ox” 
fir einen beliebigen symmetrischen Tensor M,,. Aus (59) und (60) leiten sich die 
Gleichungen fiir das verallgemeinerte elektromagnetische Feld ab 


yra My 0 É (60) 


Age o e 
He ug, | (61) 

e y—y | 
P, (mou^ U + V—y U”) — 0. (62) 


Das Liesche Partialdifferential von (53) nach den A, , € sieht folgendermassen aus 
i 1 ! 1 : a 
D [4,4] 5 = Hig |— Vv He —— pw DA, + > (moutw + 


+ Vy we) Dy | (dx). (63) 


^ 


Mit (9) führt die partielle Integration zum Ausdruck 


D LA., e S = f Vu Ee (mou ut + V—y Ur) — A, = Va Hie — 


— — oue) | wde (dx) + JA | Las Va Hus = eu) ZA œ de (da). (64) 
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Die Impuls-Energie-Tensordichte des verallgemeinerten elektroma 


l gnetischen Feldes 
ist durch die Formel i 


TW = — Yop (mo u^ uf +/—y Ur) — A, i x p. Hee — 7 (65) 
An C 


gegeben, und das Kraftobiekt des Feldes durch 
F l — e 
Vi GE An Gei Va He — gu — Du Tyè sped TE (66) 


1 
(Tay — 2 ge (O gos + de Sao — 77) 


Die Tensordichte (65) entspricht der Impuls-Energietensordichte (45) des Gravita- 
tionsfeldes. Die Erhaltungssätze für den Impuls und die Energie des Feldes haben 
die offenbar kovariante Gestalt 


Vu I. = Be (67) 


wo die Vektordichte F, sich durch die Formel 


n l — > e 
F, = irren etw) Vv Au (68) 


ausdrückt. 


Wenn wir das Liesche Differential für das dritte Glied in (53) berechnen, das 
eine Skalardichte mit dem Gewicht +1 ist, so erhalten wir 


D [Ag H af Has Hey — y (dx) = D ie Va Yop V — y Une + 


* E Hd. ie) w" de (da). (69) 
TT 


Aus der letzten Relation folgen die Identitäten 
1 ! 
A, Ur — i5 Po ZA. (70) 


die den Identitäten von Bianchi in der allgemeinen Relativitätstheorie ent- 
sprechen. 

Die kovariante Ableitung von (62) führt mit (61) und (70) zu den Bewegungsgleich- 
ungen (56). Bei Beachtung von (70) und (67) kann man auch die Bewegungsgleichungen 
ableiten. Wir sehen also, dass aus den Gleichungen des verallgemeinerten elektro- 
magnetischen Feldes die Bewegungsgleichungen folgen und dass die Impuls-Energie- 
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Erhaltungssätze für das verallgemeinerte elektromagnetische Feld identisch mit 
den Bewegungsgleichungen sind sogar dann, wenn die Feldgleichungen nicht 
erfüllt sind. 


6. Die Feldgleichungen und die bevorzugten Bezugssysteme 


Wie wir es schon gesagt haben gibt es zwei Methoden für die Lésung der Feld- 
und der Bewegungsgleichungen. Eine von ihnen beruht auf der Wahl eines bevorzugten 
Bezugssystems und der Angabe entsprechender Anfangs- und Randbedingungen. 
Im Falle des verallgemeinerten Gravitationsfeldes ist die Wahl des bevorzugten 
Bezugssystems dadurch erschwert, dass die Potentiale des Gravitationsfeldes und 
die des Trägheitsfeldes eine physikalische Grösse bilden. Nach Fock (1955) ist das 


bevorzugte Bezugssystem in der Gravitationstheorie durch die Harmonitäts- 


bedingungen 
em az 
O 
festgelegt. Die Gleichungen (43) und (71) bilden ein System von vierzehn anabhängin- 
gen Gleichungen für die vierzehn unbekannten Funktionen p* und g,,. Diese Gleichun- 


gen haben bei entsprechenden Anfangs- und Randbedingungen eine eindeutige 
Lösung. 


71) 


Im Fall des verallgemeinerten elektromagnetischen Feldes ist es leicht das be- 
vorzugte Bezugssystemm zu finden, da die Potentiale dieses Feldes zwei physikalische 
Grössen bilden. In einem Intertialbezugssystem haben die Trägheitsfeldpotentiale 
die Gestalt (51). Bei Beachtun gder Tatsache, dass in diesem Bezugssystem die kova- 
riante Ableitung gleich der gewöhnlichen ist, erhalten wir die Gleichungen des ver- 
allgemeinerten elektromagnetischen Feldes 


9, (ng w + UF) = 0, (72) 
4 
9% Heb — = o us — 0. (73) 
Diese Gleichungen bilden ein System von acht unabhängigen Gleichungen mit acht 


unbekannten Fuktionen 4, und £; bei entsprechenden Anfangs und Randbedigungen 
besitzen sie eine eindeutige Lösung. Wir sehen also, dass die Wahl von vier beliebigen 
Funktionen im Infeld-Wallace-Verfahren (1940) in der Elektrodynamik auf der 
Hinzufögung der Gleichungen (72) zu den Maxwellschen Gleichungen (73) beruht. 
Vom Standpunkt des allgemeinen Relativitätsprinzips aus stellen die Gleichungen 
(72) und (73) das Gleichungssystem eines verallgemeinerten elektromagnetischen 
Feldes im Inertialbezugssystem dar. 

Herrn Prof. W. Slebodzinski und Herrn Prof. R. S. Ingarden möchte ich auch 
hier fiir kritische Diskussionen und maucherlei wertvolle en bei meiner 


Arbeit herzlich danken. 
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EQUATIONS OF MOTION OF ROTATING BODIES 
IN GENERAL RELATIVITY THEORY 


By W. TULCZYJEW 
Institute of Theoretical Physics of Warsaw University 
(Received June 28, 1958) 


This paper is devoted to the derivation of the approximate equations of motion for 
a system of two rotating bodies treated as singularities of gravitational field. Use was made 
of the method given by Infeld (1954). Also, the equations of motion are integrated and 
a new relativistic effect is obtained, namely, motion of the equinoxes at a constant angular 
velocity. 


Introduction 


The problem of the motion of bodies in general relativity theory was first solyed 
by Einstein, Infeld, and Hoffmann (1938), who gave approximate equations of motion 
for two bodies of comparative masses treated as singularities of a gravitational field. 
The integration of these equations was performed by Robertson (1938), who showed 
that the relativistic effects reduce to the motion of the periastron. The Einstein- 
Infeld-Hoffman equations concerned spherically symmetric non-rotating bodies 
represented by singularities of the single pole type. The object of the present paper 
is to derive the equations of motion for two rotating bodies nent by en. 
of the pole-dipole type. 

The Einstein-Infeld-Hoffmann method was Bik pa improved by Infeld 
(1954). The method of Infeld is based on introducing singularities into the equations 
of motion through the intermediary of the „energy- momentum density“ tensor 


constructed from Dirac 6 functions: ; 
p D é 
Tei = Ji (à (n— & (0. ge reg 
D 


er are the coordinates of the p-th body, and d(x’) = à (x1) ô (x?) ô (x3) is the three- 
-dimensional Dirac ô function. The Greek indices run through the values 0, 1, 2, 3 
and the Roman indices, through 1, 2, 3. The equations of motion were dirived in 
Infeld's paper from the conditions of integrability of the field equations!: 


TE = 0, Héi 


1 The semicolon signifies covariant differentiation; ordinary. partial derivatives gen be denoted by 


primes. Repetition of indices indicates, as usual, sv mmation, 


(37) 
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CP !!!!! 


by the integration of these conditions over the three-dimensional neighbourhoods of 


the singularities: 


f TE de x = 9. (1.3) 


D 
Q (ST) 


From conditions (1.2) we can derive further equations 


p 
f (x — é”) Te lo = 0, 


p 

ag) 

p D 6 

li (Gi — &) (FET pdx = 0, 

er) 

A p p 

[er FE) e (an — En) TE, den = 0, us 
Qr) 


of which only the first are non-trivial, while the remaining ones are fulfilled identically. 
In the paper of Infeld (1954) the equations (1.4) were not taken into account; it turns 


p 
out, however, that the form of the coefficients t”? assumed there: 


ee egy ee 


um); pal. BIÜU I, (1.5) 


is in agreement with equations (1.4) and can be deduced from them (Tulezyjew 1957). 
The use of the ordinary Dirac 6 led to some divergences which were removed by re- 
normalization of the mass.of the bodies. In later papers of Infeld and Plebañski (1956, 


1957), there is introduced the function ó, which makes it possible to avoid renormali- 


sation. The ó function differs from the ordinary Dirac ó by the property 


iP (2) ó (xs — £)d,x —0; r= [(xs— &) (a E). (1.6) 
Qe) | 


For the ordinary à function the integrals (1.6) are not defined. Property (1.6) permits 
the definition of a “finite part“ of a singular function by means of the formula 


FO = f fe D (6 — El) dyr. | (1.7) 
AE) : 


For regular functions we obviously have 


FO = FEO, 9. (1.8) ` 


Henceforth in this paper we shall omit the A over the 6. 
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Extending the Infeld method to the case of singularities of the pole-dipole type 
we assume an ,,energy-momentum density” i 


à : A p p P D p 
T* (e — rebò, ); o = d(x — N, p=l 2, (1.9) 


and we shall seek the equations of motion in the form of (1.3) and (1.4). We shall 
use the “new approximation method: introduced by Einstein, Infeld, and Hoffman 
and applied to the investigation of the slow motions. In accordance with this method 
we expand the gravitational field into a series in powers of a small parameter A: 


3 
2 4 
. (1.10) 
3 5 i 
3 ð i 
2 4 


The expansions of the components T'9, Tor and T” will begin from the orders 2, 3, 
and 4, respectively. Hereafter we shall omit the parameter? 2 and denote the order 
of the individual terms by numbers written below. We shall remember that the veli- 
ip ? 

cities é" and the accelerations ér of the bodies are of the first and second orders, 
respectively, since differentiations with respect to time increase the order by one. 
It should also be remembered that differentiation with respect to time during the 
expansion of the field equations and the equations of motion increases the order. 
The “new approximation method” is worked out in detail in the papers of Einstein 
and Infeld. | 

Parallelly to the singularity method another method, based on the description 
of the bodies by means of the continuous energy-momentum density, was worked 
out chiefly by Fock (1939, 1955) and Papapetrou (1951). In the case of non-rotating 
spherically symmetric bodies this method gives results identical to the results of 
Infeld, but with considerably more complicated calculations. Fock also considered 
the motion of rotating bodies (Fock 1955), but did not derive equations in a form 
ready for integration. Papapetrou's method was used for rotating bodies by Haywood 
(1956). Haywood, however, neglected terms proportional to R/r, where R is the 
order of the dimensions of the bodies and r the distance between them. The equations 
he obtained do not differ in an essential way from the equations of non-rotating bodies. 

Recently there has appeared the paper of Ryabushko (1957) in which the author 
applies the method of Infeld to the motion of rotating bodies. Ryabushko, however, 
made an important error in his assumption of the energy-momentum density which 


is in contradiction to equations (1.4). 
NS 


2 In Infeld’s paper (1954) the parameter A was absorbed only at the end of the calculations. 
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2. The Equations of Motion 


In this section we shall consider the derivation of the equations of motion for 
two rotating bodies represented by a pole-dipole singularity. The “energy-momentum 
density“ of such a system has the form 


1 1 
TPO (e eet (CEE 


2 2 
+ t? à re — C3, (2.1) 


1 2 
where 9" = & and & = & denote the coordinates of the bodies. Expanding? TE B 


Tu TO ETS TO (aj in 
(en, HOA) + 240 (8) ci (9 ó — 

— re. He o 0 H2 Eg) + 
pe ta) 5, enge re ds + 


+ similar terms for the second body (2.2) 
we obtain from (1.3) and (1.4) the equations 


— 1 


1 CA aie 1 EEE, 
RT Oo rn. Ce ek 


TAN II Mere i 
se oM PAM He fed CS 2. 3a) 
1 1 1 1 EL 1 ta 
AA ab WÈL ee TA kan Aa cet, Me pa AR dual, (2.3b) 
1 1 1 1 : 
Gen 7 — 1560 7 “= pe cU per = 0. l (2.3c) 


We obtain similar equations for the second ‘body. The last of equations (2.3) 


can be solved exactly, since they do not contain the grawitational field. Let us first 
1 2 


assume that 1799 and 70 vanish. For o = 0 equations (2.3c) now give 


1 1 | 
Poses OMe: (2.4) 
Wie introduce the notation | 
i 1 1 1 i | 
vo er tor =S (2.5) : 


è (a) = 3 8? Leon, „ + go, 8 — Epy, o) denotes the Christoffel symbol of the second kind. 
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From equations (2.4) it follows that 


1 ó j 3 
TALAN pon LAJ 
Dar oo ei: (2.6) 


From equations (2.3c) for & = t we obtain 
1 11 1 
5 S* + 3 Styr — pres — tr = 0. (2.7) 


Adding and subtracting from these equations the respective equations resulting from 
a cyclic change of the indices r, s, t we finally obtain 


1 je A, lae 
prope ES (Ss + Sra). (2.8) 
We introduce another symbol 
nor 
ni == 188, (2.9) 


Equations (2.3a, b) can be solved only after the field equations have been solved. 
We take the field equations in the form 


RP — — gg p*? (2.10) 


where R7? is the Ricci tensor density and 


pe = a = gaß gs ed. | (2.11) 


We shall need the quantities hog, Amn» Pom aud hoo for the equations of motion in 
2 2 4 


the first post-Newtonian approximation. The expansions of the field equations given 
in the Appendix lead to the approximate equations* 


11 2 2° 
Ahoy = 87 (m ô + m ó); (2.12a) 
2 2 2 
Ab, = Om Aho (2.12b) 
2 2 i 
1 de ew 2 2 
Ah, y = — 167 EE Sm ó, , + töm ó _ $m a] , (2.126) 
3 3 3 2 3 
3 1 ek THUR ii im. 7 2 2 Y 
ET P rr — STs as ir OM 
dn Bri ti )ô é n Mu d 
DE 
SS 650, r] SE hoo: 00 Se Roos Foo: s 3 
3 2 2 2 artes 


3 1 
ca eh EN (2.124) 
— hers Roos rs + 2 500200 t 2 5 ARA | 


4 In these editions certain conditions were imposed on the coordinate system (see the Appendix). 


We shall return to the problem of the system in the next section. 
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here A denotes the Laplacian operator. Equations (2.12 a, b, c) have the solutions 


1 2 2 Pp p 2 
hoo = — 2m/r— 2 Ee r= [(x° — E (2° — EN (2.13a) 
2 
hon = Onn l op: (2. Läb) 
2 2 
för al 222 1 1 2 2 
hom = 4t""[r +400" jr — 28" (l/r), — 25” (Afr)... (2.13c) 
3 3 3 3 3 


Employing in equation (2.12d) the results obtained thus far and making the transfor- 


mation 
va 
hoohoo,ss + bag Roo,s = 2 (Loo) , 
1-1 2 2 
hoo hooss = 8% hoo (m ó +m ò) 
2. 2 
1 2 
mm è 
= — 161 (ö d ò) IG A-n (2.14) 
we obtain 


1 1 1 1 8 l 2 2 2 2 Sëch 
A hog = 82 [(m + t") ő — Sr ő, + (m +07) ó — SA öd, — 
4 4 4 3 4 4 


Ex (m min 6 ＋ 4- hong +5 Die (2.15) 
from which we have 
1 1 1 2 2 2 1 1 N 
hop = — 2 (m Her — 2 (m H fr 2 (Une + - (216) 


2 ok 2 172 1 i 242 
+ 25" e Un, SCH no) (/ + +n) +2 (mir + + minè — m (oo xg (7),00- 


We may now proceed to solve equations (2.3 a, b) and eliminate the gravitational 
field from it. We write the approximate equations in the successive orders of the 
approximation. For & = 0 only the odd terms occur, for & = m only the even ones. 
We split equations (2.3b) for x = m into parts that are symmetric and antisymmetric 
with respect to the indices r, m. We begin with the third order 


= 0, 
Tg (2.17) 


tor = mi. EN (2.18) 
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In the fourth order we have 


nim Ale) 13 nm (>) 
mn" + m =mn"—mm|— ; 2.19 
1 2 a; 2 d E NI Fo Sa 
1 19 
en m in , (2.20) 
4 2 
Zen, „0 = 0. (2.21) 
3281 


The expansions of the Christoffel symbols are given in the Appendix. Equations (2.18) 
and (2.20) are or an auxiliary nature. Equations (2.17), (2.19) and (2.21), together 
with similar equations for the second body constitute the Newtonian equations of 


1 2 1 2 
motion if m and m are interpreted as the mass of the bodies and S“ and S" as the 
2 2 3 3 


internal angular momentum connected with rotation. We now pass on to the fifth 
order. 


1 1 01 1 0 Wat 9 1522/1 
ng tepore mm he ef (7) r= 
3 > 
TELE 18872 
Im i Nn m mir , (2.22) 
i 2 2 2129.9 
Thus 
: 2 
m EIN n” +m mir (2.23) 


In (2.22) we made use of the Newtonian equations of motion and of the transfor- 


aer 1 . 1 
4 + (=) ll ge Gi 


In (2.23) we omitted the integration constant, which can always be added to the 
Newtonian mass. Further equations of the fifth order are . 


mation 


1 241 il 
1 ES 1 DIETE TEL | 2.95 
re), E 


We now consider the equations of the sixth order, that is, the post-Newtonian 
equations of motion. From equations (2.3a) we take the terms of the sixth order: 


ton,, ＋ mim + mim + 2m nf (23 + 
5 2 4 4 2 2 = 


133 E pe LAN: Ge 
+ man S" Veer SMA (e Ze 
f 2 2 3 3 3 2 
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These terms should be regarded as corrections to the Newtonian equations of motion. 
The requirement that these terms vanish separately would lead to a contradiction 
with the Newtonian equations. Equations (2.3a) for & = m can be satisfied in each 
order separately only if the coordinates of the: bodies , and D are expanded into 


an approximation series. For x = 0 equation (2.3a) can be satisfied in each order, 
1 2 
since masses m and m were expanded into orders. After a suitable transformation 


of terms (2.26) and addition of them to the Newtonian terms we obtain the post- 
-Newtonian equations of motion?: 


1988 ik Sat ki A Forge SE 1 
mn" —mm B =m 1 7 ＋ 18 fs — An! b —4＋ — 
sm / 


1 

i . . . . . E . 1 

are B + [Ans (£n — n”) + 3 n ès — 46” Cs] 5 Ss 
Sa ‚m 


1 TA E : 12 : : 1 
+g Orem mer — 28) + S emis 


aui d 19 AG. a 1 
+ [2 m Sr (ns — £3) +2m Sr (n° — CH B (2.27) 


nnt 
In these equations we kept only the terms proportional to R/r, while omitting the 
125] : 
term GE AC ( Jp proportional to (R/r)? as being considerably smaller than the 
"n fí n" 


r 


others. The equations of motion for the second body can be obtained by changing 
1 REI Rd 5 

n, C, m, m, S", S" to &, «y, m, m, S", S" respectively. For simplicity we have 

omitted the indices denoting the order of the approximation of the quantities 

é ' 


P 
m and 5”, 
D H H D | . 8 á P» D 
We have not yet justified the interpretation of the quantities m and S"? occurring 
in the equations of motion as mass and angular momentum and the assumptions 
1 2 
. : (2.28) 


which we used in deriving the equations. The singularities represent the gravitational ` 
field outside the body, that is, in the region in which the equations for empty io 
are satisfied: 


R | (2.29) 


5 We will not give the precise transformations leading to equations (2.27). These transformations 
are rather lengthy and to a large degree are repetitions of the transformations contained in the paper 
of Infeld (1954) on which we continually base ourselves. 
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Inside the body we have 


ki 1 
Re = ön (ou n. su, or on) (2.30) 


where O9? is the continuous energy-momentum density tensor. The equality of the 
gravitational fields obtained from equations (2.30) and (2.10) outside the bodies is 
the basis of the interpretation of the coefficients tab and pe. Let us take the equations 
in the lowest approximation: 


2 2 
Ah = Ba, (2.32) 


The equality of the fields hog is obtained by requiring that the multipole moments 
2 


of T' and 099 be equal, that is, 


2 2 
1 1 o 
00 m = f G00 d, 100 m = f Ody x (2.33) 
: oe) 2 ” od» 2 
and 
1 2 
7700 — f (x? — 17) Od. x, 1700 = f (a? — £7) G00 % x. (2.34) 
2 Anm) 2 AINE pure) 2 | 
The higher moments do not interest us because of the approximation in R/r which 
is used. O is the Newtonian mass density; we thus see that íz and = are really the 
total en masses of the bodies and goo and goo are the static moments of the 


2 2 
bodies with respect to the points with the coordinates 7" and &“. Conditions (2.28) 


in the Newtonian approximation therefore denote the fact that points y” and & are 
the centres of mass of the bodies. We shall examine further the equations 


Aho, = 16 OU, (2.35) 
3 3 
Aho, = — lén T”. (2.36) 
3 3 
Here too we should require that the multipole moments be equal: 
1 Sek 
o = f Ordaz, = dh (2.37) 
3 Ans) 3 3 ags) 3 
and 
1 2 
ep — G 1 Erd *. 4 — = Mr dg x. 239) 
3 am 3 3 ay 3 | 
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1 2 
In this way 8“ and S" obtain the interpretation of internal angular momentum, 


or of angular momentum computed with respect to the centres of mass 


1 
Se =f [67 — 77’) I — (x! — n’) sg? de x. 
Ane) 
2 e 
$8" [elte L MOU — (* — DON dya. (2.39) 
Qs) 2 3. 
We have given the interpretation of conditions (2.28) in the Newtonian approximation. 


In deriving the equations of motion we used these conditions also in the post-Newto- 
nian approximation 


1 1 2 2 
goo 1 pr00 = 0, goo 42700 — 0. (2.40) 
2 4 2 4 


If we were to replace conditions (2.35) by the Newtonian condition 


1 2 
ro — 0. 70 — 0 (2.41) 
2 2 


we would obtain on the right-hand side of the equations of motion additional terms 
which we give without the calculations: 


1 21 Lio entre 
— Sm. % + (m S —m S") | — , (2.42) 
AMES m 

1 2 1 2 
S” and S” denote £709 and 499. Equations (2.27) thus refer to points displaced with 
4 4 
respect to the Newtonian centres of mass and differ from the equations of motion 

1 2 

of the centres of mass by the terms (2.42). The quantities S" and S” are in general 


dependent on the relative positions and velocities of the bodies; a precise interpreta- 
tion of these quantities is difficult, since the equations for hy, are rather complicated. 
4 


In the next section we shall show, however, that terms (2.42) can be made to vanish 
with a suitable transformation of the coordinate system without affecting the rela- 
tivistic secular correction to the Newtonian equation. The interpretation of the quan- 
tity 5" thus becomes non-essential. 

After integration of the equations of motion it is convenient to use the Lagrangian. 
The Lagrangian of the system of non-rotating bodies is known to be (Infeld 1957) 


1 2 
T Jw 
2 tpm gå Cy he dez (CHAN 
1 2 1 
l mm i — 
S Se scs en) 
1 2 
lmm 


2 gc (-Uwm-ob. €43) 
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The full equations, with the rotation taken into account are obtained, as readily ve- 
rified, by adding to (2.43) further terms | 


1 2 
d ce é is *s r 7 d s 2 = 7 
% meer (2.44) 


3. Integration of the Equations of Motion 


The derivation of the equations of motion in the preceding section holds in 
a coordinate system in which the gravitational field satisfies the conditions 


Sap rz Nag? Mo = l, 7], — Lori o No = 0 (3.1) 
and 
l 1 
ling r h 55 FT 
l > zs 2 zou 2 i m 0 
l 1 
or, EL 2 hooo — 2 frr, = 0. (3.2) 


Condition (3.1) denotes that in the zero approximation the most natural coordinate 
system for a flat space-time is taken, namely the Galilean system. On the other hand, 
condition (3.2) has no physical significance, and its adoption can be justified only 
by calculational reasons. Passing to more general coordinate systems we. shall use 
a transformation preserving (3.1). If we require additionally that the transformation 
not violate the assumptions on the parities of the orders of the expansions, then 
we obtain® 


x — x" +a Lat 
2 4 

PLC YE TEE (3.3) 
3 


where a” are arbitrary functions x". Inserting 


ér u Cz se a’; = f a’ ó de x, (3.4) 
i 2 2 


into the eguations of motion we obtain on the right-hand side additional terms of 
the sixth order ( 

SA i Ag ie tte 

— m dio + m m (ar — a") ( (3.5) 

Zë 2 CA ge gm 
cac coe Coa au Re —T—-—ꝓ6 EE ae 
* We neglect here space translations and rotations; these transformations will play an important 

role later on in this section. 
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Functions a” and a9, on the other hand, do not affect the equations of motion in 
4 4 
this approximation. Let us now consider in which of the systems available by means 


of the transformation 


N (3.6) 
2 


the post-Newtonian effects observed by a distant observer should be calculated. The 
results of these calculations would obviously also be the same in all systems, but we 
shall attain our purpose more quickly by imposing the following conditions on the 
system: 


(a) the gravitational field depends on time only through the coordinates of the 
bodies 


(b) the gravitational field takes the Galilean values at infinity. 


These conditions can, for sure, be satisfied; for example, our initial system 
satisfies them. Let us suppose that we have integrated our equations of motion in 
a system satisfying (a) and (b) and that we obtained, as a result, periodic motion. 
We maintain that the motion observed by a distant observer by means of light signals 
is also periodic and with the same period. Indeed, if the bodies return periodically 
to the same positions with respect to thé system and simultaneously by condition (a) 
the metric of the world attains the same value periodically, then the results of observa- 
tions will also be periodic functions of time, since the signals sent out by the bodies 
will be emitted periodically from the same points and in the same conditions. Condi- 
tion (b) ensures the agreement between the period of the observation and the period 
of the motion with respect to the system; from which it follows that the time coordinate 
coincides with eigentime of the observer resting at infinity. Systems satisfying condi- 
tions (a) and (b) are, as may be seen, very convenient to investigate periodic motion, 
particularly if the period of the motion is the only quantity in which we are interested. 
By using such systems we avoid lengthy analysis of the light signals. The motion of 
the bodies observed over a short time is described, to a sufficient accuracy, by the 
Newtonian equations of motion, the interpretation of which does not entail any difficul- 
ties as regards linearity and constancy of the speed of the light signals in the Newto- 
nian approximation. Relativistic effects appear with slow periodic? changes of the 
parameters of the classical motion. The relativistic effects can be calculated by the 
step-by-step approximation method. In various systems we would, in general, obtain 
different results; but if we carried out the integration in a system satisfying conditions 
(a) and (b), then the calculated periods of the secular changes coincide with observed 
periods. Our initial system is not the only system with these properties, since transfor- 


mation (3.6) preserves conditions (a) and (b) if only the functions a” vanish at infinity 
2 
— T—HVT:Pt Car .. 


Non: periodic changes may appear in connection with gravitational radiation. The radiation doe 
not affect, however, the equations of motion in the first post-Newtonian approximation. | 
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and do not depend explicitly on time. In particular, setting 
1 1 2 
ma =—S, ma =— S" (3.7) 

2 2 


2 


~ 


1 2 
we can absorb terms (2.42) without violating conditions (a) and (b), since S" and S" 


can depend on time only through the positions of the bodies. The discuccion given 
here is in agreement with the results of Haywood (1952), who transformed the equations 


of motion of non-rotating bodies from the De Donder system to a system satisfying 
the condition 


gg = 0 (3.8) 
and showed that the period of the motion of the periastron in both systems is the 
same despite the fact that the shape of the equations is changed in the transformation. 
Both systems satisfied conditions (a) and (b). 

Let us pass on to the integration of equations (2.27). We use for this the Lagrangian 
derived in the preceding section. Introducing a number of new notations 


a F , 
1 2 
A = ag (mm eh, 
1 2 
M = m n, 
2 1 
Om == mm 
1 2 
mm 
AE Sv 


1 2 2 
Mòs" = (1 + m/m) S" — (1 + m/m) S”, 


1 12 2 
M2[k, l] = km? + Imm + Am, (3.9) 
we rewrite the Lagrangian in three-dimensional vector notation: : 
il 1 1 Ai e e DESIT 
„% [[[ MKR 
1 uM 1 pòm =. 1 uM san 1 pdm „7 75 
ay = e 2 2 po g r? (n Dar ALTE 


VVV E IA 
2 = SITE BEE 


Ke LEE 3.10 
e e s 
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The vectors R and 7 have the components X” = xm p R, D = 7, the vectors has 
the components d = 4 e" s" A 3 NL =. F r / / eke 
= 1) similarly with the vector és. 

Making use of the invariance of the Lagrangian with respect to certain transfor- 
mation and using Noether’s Theorem (1918) we derive a number of conservation 


laws’. We begin with the displacement 


N Ret dR (3.11) 


where OR is an arbitrary constant vector. Noether’s Theorem gives in this case three 
momentum integrals: 


25 1 H 1 c "o 1 pòm > uM > 
i My? — uv ? 
N MV V Lam ae M 
IA 7 1 
wa 3 Ging * . (8.12) 
r r 2 


From these integrals we derive the conservation law for the centre of inertia 


(mature HM) K (a fiat are 


DENIS Pe don A 
—p (9 x 63) = Pt +k. (3.13) 
In the system connected with the centre of inertia (P = 0, K = O) we have 
> 106m [uM s 
wR- 5 (45 ; lee) ap x 63) (3.14) 


Restricting ourselves to such systems we can discard in the Lagrangian all terms 
containing the vector V as small at least to the eight order. Neglecting the factor u we 
obtain the Lagrangian of relative motion 


ana 3 


2 l, 11 -7 Ze vye re AEN (3.15) 
From the invariance of this Lagrangian with respect to displacements in time 
t t + ôt : (3.16) 
there follows the energy integral 
1 M 3 1 mee 1 M 
2 „ l. % 4-3. [39] — LESEN Ze pcr Dig o5 
(3.17) 


8 The use of Noether's Theorem to derive the conservation laws in Newtonian mechanics can be 
found in the paper of Hill (1951). 
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Invariance with respect to rotations 

er, 

„5 (3.18) 
gives 


> ES —-— mm 
beer Side = 2-5 $ xi J, (3.19) 


where 


E $ l NE 
i= 6x9 (14 2 U. . .) +2; x (r X 8) (20) 


The scalar product of (3.19) and i gives 
J = | J| — const. (3.21) 


Using the energy and angular momentum integrals we shall find the path of relative 
motion in typical special cases. Let us first assume that the motion in the Newtonian 
approximation takes place in a plane perpendicular to the vector 5. From (3.19) and 
(3.20) we now obtain 


J—Gxö) (2 nep x 2; weg (32) 


r 


In this case the motion is plane also in the post-Newtonian approximation. Introducing 
polar coordinates in the plane of the motion we obtain 


M M. 

J=rġ (ien nes n) SEE Qo 

1 MES 1 M LM! Me «es 
EE tglh — 1] vé +3 B Sp =: 5 ==: lu (r? oli, 
(3.24) 

Determining from this v? and r?g and inserting them into the identity ` 

CAW 2 12 (2 9 ; (3.25) 

dp +u? = vè (r?° p) : 


where u = l/r, we obtain 


En Im ra[ag +0,18 + (2M + [8, 10] em — 8-5 EM) NE 
p 


«dte 9] M? — a = we) ul + [0, 1] Mu? (3.26) 
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and by differentiation 


d? : E A EE 
doi +u= JAM 8 EM — 47 EM - (l. we — 83 me) uf + 
3 1 M 5 s 

y (0, 1 Mut = — perd e (u. 5] —4 ai +m (16.9 8 a ul | 
+ [0,1] Mu? (3,27) 


where E = —M/2a, J? = Mp. In order to obtain the periodic solution? we put (3.27) 
into the form 


dën 1 1[M 
dg? 8 t à mr e (i es 13) SN 


m (169-85) : d ent [0, 1] Mu? (3.28) 


where a is a quantity of the second order undetermined for the time being. In the 
lowest approximation we obtain the solution 


* = (1 + € cos (1 — a)g). (3.29) 


Inserting (3.29) into the right-hand side of equation (3.28) we obtain the equation 
of the next approximation. This equation will have a periodic solution only in the 
case the resonance term 


a Q 9] NE + [0, a) 27 cos (1 — a) @ (3.30) 
vanishes. From this we obtain the condition for a: 


M s 
5 6430 (3.31) 


Further solving of equation (3.28) is without purpose, since & is the only quantity 
of interest to us in this problem. The solution (3.29) represents a rotating ellipse 
with the periastron advance by 


M 4 s * 
À © = 2na = 6x — (1 — — 3.3: 
or 3 | 3 A (3.32) 
for the period of classical motion. 
——— I Sa c ueteri Were a ĩ·¹.w... ice EE e 


° The method used here is presented in the book by Krylov (1950). 
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Let us next examine another case. We assume that the Newtonian motion takes 
place along a circular path r = a. Equation (3.19) now gives 


ip Ane E (3.33) 


The vector J which is perpendicular to the plane of motion now precesses about 


> 


a fixed vector s with an angular velocity 


Q2-— s. | (3.34) 


This velocity is simultaneously the angular velocity of the line of equinoxes, that is, 


the line of intersection of the plane of motion and the plane perpendicular to s. The 
results given here are in agreement with the results of Lense and Thirring (1918) 
for test particles moving in the field of a rotating body. 

The author would like to thank Professor. L. Infeld for suggesting the problem 
and for valuable advice. | 


Appendix 


The expansion of equations (2.10) into an approximation series 
On the left-hand side of equations (2.10) is the Ricci density tensor: 


RA = g gg R, g= ga (A. ]) 


We expand the individual factors 


x i j 
V—g = U 8 1 e ee à (4.2) 
We obtain the expansion of g% from the relation 
fy Ve 5 ' (A.3) 
Ben „„ D atv : $ 


It may readily be verified that 


pos = 1 vo} hoo am — 
2 i 


535 (A.4) 


For the expansion of Rag 


Rag = (ha — fashe. À lan) {oa} — tap) Go) (A. 
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we shall need the expansion of the Christoffel symbol 


k 1 1 1 
ee en EE ling ee 
loo 9 100% t3 100, loko + 9 ler 1200, SS 


1 1 1 1 
Ro = — — FES ae ey AG = Gy LoS Ti 
00 9 QUE 2 Tons ^ Ross 2 hss,09 2 2 $3 hoo, 
7... ee 
2 à rs 1 00,rs 4 a 00,s a SÈ 4 a 00,5 quu HE) 
1 1 li 
Ron = 2 Pons TE 9 Den SCH 2 Debt XE 2 heni SC > 


Collecting our results we have 


1 1 1 1 
R00— —— — — IS ur pics M 
2 Dat, d 2 Lob. as hos,os 2 2 pen 2 hrs,s hoor 
1 1 1 i 
FFY 2 Ars hoo,rs + 4 hoo,s hoo,s + 4 hoo,s Le, + 
* 2 2 2 2 


: 1 1 
Rm — — — — — ern — = 1 
9 Dono an 9 Pons ns är 9 ham DE Brenn 2 hoosmn 1 


(A.6) 


(A.7) 


. ASÈ) 
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The expansion of the right-hand side of equations (2.10) does not entail any difficulty 
We obtain i 


1 1 
po. = GMD AL m T99 L Trr 
2 2 I 9) 6 i a ) + 
po» — Tr) + f 
3 
1 
pn — — Omn Zu + (A.9) 


Using in the equations 


2 2 

Rm mte AE. 8 Pr 

2 2 

R” — 8 Pon, (A. lo) 
3 3 

R” = 8 po 

4 4 


the expansions (A.8) and (A.9), the conditions for the system 


hor, — 5 hoo,o — z lao = 0 (A.11) 
3 KS 2 1 


and the earlier results of section 2, we obtain equations (2.12). Contitions (A.11) are 
obtained from the expansion of the De Donder condition 


(V— g8), p = 0. (4.12) 
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The canonical transformation which excludes from hamiltonian the interaction term 


He is analysed. It is proved, that the transformation, applied in approximation of 
first commutators with generating function of transformation, leads to changed subsidiary 
conditions and new hamiltonian expressed in the Bohm — Pines representation. It is 
pointed out, that because of the approximation degree used in applying the transformation, 
the set of supplementary conditions and the new hamiltonian must be considered as the 
equations which are rather postulated than correctly derived. There are also obtained 
generalized dispersion relations. 


1. Introduction 


In the part I of this study there was formulated a generalized form of the Bohm- 
-Pines theory for electrons moving in nearly periodic field of crystal lattice. It has 
been shown, that several expressions treated by Bohm-Pines as small quantities 
become very. large with aproach tpo the form of theory proposed by these authors. 
It has been also pointed out that the separation of the wave function in product of 
a form depending on feld variables and a second one describingmotion of electrons 
interacting by means of Coulomb short — range forces, is not exact, except in 
a case considered by Kanazawa (1957). 

In the present paper we propose to derive a new representation of the hamiltonian, 
called here BP representation. It can be obtained by the exclusion from the hamilto- 
nian (2.1) of the term Hie by the help of suitable canonical transformation. In 
consequence we obtain an energy operator which does not contain this interaction 


. member. 

Since the transformation is applied not exactly, but only up to the first commutators 
Vith generating function of transformation, there arises a situation that both the new 
` hamiltonian and new subsidiary conditions have to be considered as being rather 
postulated than correctly derived. 

It seems to be quite fruitless to seek after a more correct form of the Bohm-Pines 
theory than that obtained by these authors, because the equations of such an exact- 


(57) 
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theory, including a set of subsidiary conditions, would be not in the least easier to 
solve than the equations of a straightforward theory of electrons, founded on many 
electron Schrödinger equation. 


2. The hamiltonian of electrons in crystals in BP representation 


We shall start from the hamiltonian and a set of supplementary conditions derived 
in part I of this paper. We have 


A a A A A A 1755 
H= Hpart an Hs Ee Hoon * Hinter * U— 2» 4 85 x 40-rx) Chwa) — 
2nZ (1 + Z) Ne? 1 
Xu quet PARTIE" (D 
d ASA. 
with 
no -Å D32m + Y B2/2M, (2.22) 
La 
ho 
= X [^ (az az +a,ai) — 
i AA : | 
Ah (wi — wp) 
eres 2 (aia; + a, ai — 14 AZA and a)» (2.2c) 
Lac, 
a eh V RE BE 
Hinter = = H D? ha" À [ex (B; — h À/2) eti az + az e77 "i ez (Pa —kA[2)] — 
er 
Zeh ` 
€ ze A Ka n e. ch ^ E 
DÈ D ha A len- (Ba — KAJ) re ay tate rat, (Pe—hd]2)], (2.20) 
a=1 J, 181. 
eee SÈ 25 6 E cure 
== Ee ess dy — Ao Li g Ti 
V ; e ben 0)” a, — a" 1 o) = x e + 
3,6XA, 
A+5#0 
PO ak c 
EY an] * 
a=1 
aa, kapa ESVE Os AA, (2.3) 


All the designations are the same as in part I. 
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We know already that (¥*, U Y) = 0, therefore we treat U as a small quantity 
and neglect it on base of random phase approximation (Bohm and Pines 1953b). 
We shall now try to transform the hamiltonian (2.1) using a canonical transforma- 


tion generated by a function S given by 


$ = Sa + Se (2.4) 
where 
S EN ear (D; — hål : 
Kë = — i E ha A | = (D; EM er; he 
eier; wi A P;/m + h1?/2m 


— 03 gi T; SR (P: Ta hA/2) | ; (2.5a) 
wa — A+ m + m 
Zhe V & (Pe — h 3/2) e 
Sin = i SA? 35 vil e — - E ete ag — 
u a=1 J. 1<4, Wa — À* BelM + hi?/2M 
Er DANTE | (2.5b) 
03 — A+ PM +hi2/2M 
with 
27 á 
De 5 | 
All the operators concerned transform according to the scheme 
a = i y 1 Aa 
FVV E (*) op Lá, S], . (2.6) 


p=0 


Here [A, S], signifies p — th commutator A with S and [4, S], = A. T 
The series (2.6) converges rapidly if the parameters ey and e, are small compare 
to unity. These parameters are defined by the relations 


ue N. (2.73) 
DA E 


where the averaging is Vet through with the help of ground state eigenfunction 


i-Di istributi ction near 
of the hamiltonian (2.1) or, less exactly, for Fermi-Dirac distribution functi 


absolute zero and over sall 2 directions in half a sphere A « å, 
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Whether the nuclei of the ion — cores are besons or fermions does not play any 
róle here, since the velocities of the ions are much smaller than these of electrons 
and therefore €, « ei. 

After averaging for Fermi distribution we obtain 


de? 1 40 
| Mes: EO 2.8 
a IS a (2.8) 
SEM 
A hassen 2mus i2 
with ag = — >» PI (4nr?9/3) 1, Ay = £ where wo is the electron 
me 


chemical potential. 
. . m a], 
If we assume that the ion — cores are fermions, then e, = M Z ^&. However, 


for the present, we do not take this fact into account and proceed to estimate the per- 

turbation terms for both parameters €, and e, simultaneously. It is obvious that in 

this case we omit some electron quantities which are much larger than the ion expres- 

sions left. We neglect therefore small ion terms towards the end of our analysis. 
The following relations are easy to verify. If we assume 


AN 
4 1 1 1 à 1 * a 
H. = =— ) Bi +5y ) Pa + E. d ho, Laag: + anan), (2.9) 


i=l a-1 di 
we obtain 
Da RE H 
D Ha, S ] — Ib > (2.10a) 
therefore 
L n & A i n—1 A x 
h [H,, S]. ES e es Shes , (2.10b) i 


what finally leads to 


ei (Ha SS Hinter) ei. = Ha + Hinter + ys E (i) [Ha + Hinter > S]. 
; n=1 ` à 


x 1 A a S a A 
PE 6 ) H- (n = 1) eim | Hate. Sb = He T 3k [Ho Si + ... (2-11) 
We shall prove that (2.10a) is true. For this purpose we notice that 


[22, e I = 26, hÀ- (D, — h42) Än | (2.123) 
Lë. e] = — 205 & e. A. SZ (oc aab) 


Modified Bohm and Pines Collective Description 


61 


[Dés € a = p 20,5 hå- (p, — hide? "a, 
ES ere] = 203 he" ta J. (P,—hÀ/2) : 
[af a,,a,] = la, , a] —— On Qd. 


[ay a. aj] = [a à, aj] = 9 


A ST. 5 ap D; — hA/2 
a [A 5] a = y hit 4| SA (P, / à 
h =l w — A Bj[m +ha /2m 


D he Bim +h 42/15 


> t A 2] E ; sè 
VCC 


Wi — A: P;/M-+h A?[2M 


pen 1" (Ps — hA/2) E 
ai Li, eng] c — À: BAM +h 222M 


(2.120) 
(2.12d) 
(2.12e) 
(2.12f) 


Io 1. a, — 


[ac ag + do as aa] e "a — 


eh V HE (P. — h 4/2) z * LS 
TRAIN h ES n n El (ke -F Ag Ag, aa] e 2 
TOUS X ^» 23 De Bim + h 42/2m i 
i=1 A ASA, 
9 A. 
- ig ] “ther; ea (Pi KAY = 
alg da. te de, 0a] € w, — A+ Bim +h A?/2m 
Zee en (Pe — h À/2) 
— 3M Er DM CAEN 
2 el Lacs WI : Pa] + 
0,04, 


las as + do A, aj] e “A Ta 


EN E. — ^ ij] 
Rx T 
a i=1 2, AZA, w, — À: Pilm + h A2/2m 
MTM 

gin wi — A: Pm + h A2/2m 


pr SLR SHE 
wi — Å: B.M 4- Mu 


el^ rí a + 
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FFV 
Se wa — A: PalM + h A?/2M 
| 4 [A+ (b, — h )) J | 
„ gig DEREN ee CUIR EL A 
xa aj OA w, — À: BAM +h 2/2M 
EC AY n LOL SE 
Lp emt? ër w, — A: Pm + h A?|2m 
ale a; re eb r; Ar (D; == fj) A[2) 
Wi — À: Dim + h 2?/2m 
n 15 san 2 ^ 
a 2L ps haor | le e 
a=1 145A, wa — A: BAM + h A*|2M 
A: (Pa — h 112) | eh V EE 
-id.t, NUS US) DEUS 
he EE BAM +h A2 m M » AM 


c > 
i=1 2, ASA, 


HN) a4 I di e. 1 . . —&3p) + À yo > ha À {ex (D,— 


a=1 A, ASA, 


E DÉI FINA a, -— a; eit re & $ (Pa == h Àj2)} E es Hinter e 


We compute now 


3 Himes 


WEE Bin = purs | ne 


i ^ ^ 
25 Hinter- S] explicitly and get 


A | 2 = 
5] = » ho (a4 a- + a5 a-;) P Ds E z 
i=1 


A, ASA, 


vi = * ne? 
(o +Ï: PIM)? —h? 19 à RE | mV 


LE 


Wa nZ? e? N | 
(wi — A: Bim)? —h? Sc + MV > = 


wi 


Sarees cli 
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me L1 D; hA SDN 
-—— X LE ei (Pi + háj) ea" (Bj — ni) | 
iZ XA, wa (wa — A+ P;/m + m) 


L5 (Bi +h AY er (P; —h åf) us d 
wa (03 — A+ Pm +h 4% 2m) Ars 


LE X fane (Ba + h 2/2) ea (Pg — h A/) 


ane a+b 1121, wa (4 — A: P;|M + h22/2M) 
£i (Pa + hAl2)er - (Ds — KAB). Br «| oe 
wa (wa — A+ BM + h 22/2M) (2.13) 
2 > uten Ker?) RP 
d JET ISA (o, 4: Dim)? — h? in 
8 G: ba)? À 
VA g—1 A, ASA, (wa — À: 5. M)? — — h? 434 M? 
2 Zei X^ V 
tmr D E dr A AKA oe 


c 


= d 1 : 
2 LU B e, mm afl 
/ wi — (A: Pm + hA?/2m)? ` 9$ — (å Pa| M — hå? /2M) 
In the above expression certain terms are omitted on the basis of random phase 


approximation and several others on ground of their proportionality to [m p bh 
mwa 


> À 3 
or a and higher powers. 
M Wa 
When we suppose that the. dispersion relations have the form 


Ame? V 1 
mi pce YE EE E 
my £4 (co, + A: Pim)? — h? 14/4m 


N 
SE pacc Gor sved SVA (2,14) 
MV £4 (wy À: Bal M)? — RAM 
then the expression (2.13) becomes 
i ee a Wi — Vè yn . EN 
op Minter, S] = Kane Lag — aa — 434 3) — 


> 
Aish, 
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me? EZ Er (B; +h Â/2) ez: (P; — h A/2) 
m?V Zu i Wa (wi — À* ÉiJm + h A2/2m) 


q Ex (Pi + h 22) aa: (Bj — h 72) en] 3 
w, (Wi — A+ Dien + A A2/2m) 


= * [e rag Sa (Pa EAD) E. (Bp — h 1/2) 
NE reine w, (01 T. aM + ri 


+ 


en (Pa + hA/2) & (Ds — h À/2) Jis 1 
wa (o; — Å- P/M + h A2/2M) 


EE Gär 


mV e 1, (o, — Ñ Bim)? — h224/4m2 — 


2x Z?e? V (& PJ? 
=e — — ch 
wr À, 2 (.J. B. Mae — h22 JAM? 


N 
SS e 
i Aal, 


1 1 
—hi dL o NI 
a) E — (Bim , wi— (PM —h 2M) ds d 


n 
ne? 


m2V 


The terms — 


. 

1 XA ( — J. Pm)? — k? 34m? A 

haze À Gän 

— Z, Ds 2 rs de , Can be casily transformed 
MeV ee (w — J. 5./ N) 2 — N 

into the form 


. 


Mrs 92 Le 50 då = T 
i=1 


n 


n’ 
Ka V 3 
? Gett 


D i 
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n' LT ` No L 7 AB V a 
a KE py sop ee e y musee 1 pans 
3n | A p / m;n (223 | di 6772 m (2.16a) 


i=1 TE A, les 
and quite similarly 
N H bé 
= Ser) 3 3 (ea . Pa)? T Te 7 m M PE/2M. 
dck? a=1 À, 1<4, (wa — 1-P,/M)? — h2 AM? 3n M a 


(2.16b) 


It is still necessary to transform the dispersion relations (2.14). If we assume 


N 
that the ion — cores are nearly at rest we get >) B, = 0. On the other hand we guess 
a=1 


ñ 
on ground of the momentum conservation law 2) f; ~ 0. Hence 
i=1 


Are? d 1 


55 AE PER > A > A Se 
mV ém wi {1 — [F 2}; Bimo; — (F: Bmw)? +(kA2mwa)2]3 


AnZ*g? À 1 2 
MV 2 oi 21· P. M — (I. / Mo;)? + (h 22/2Mo;)?]) 


-= y 1/03 [1 F 2 (A+ Bi[mo;) = N P. + (h 22/2mo)? + ...] + 
mV | A 


N 
2,2 A St 
an Ke 8 1/03 [1 F2(X-P,/Mo,) +3 (4 : Pe Mo)? + (h4/2Mo;)? + . J. 
vå (2.17a) 
and 
3 AC EVE Ze Raye árta 
ar G. 50 As L ÄB) + 
n: m? (1 +27) i=l a=1 
Le A aide: Ts : DA , 7b 
+ mr YO + 2 (4 Bi)? + 4783 (2.17b) 
where wp = ZEN denotes the Langmuir frequency cf electron 


mp ME ing i f the crystal 
oscillations corrected by taking into account the slowly moving Jong pe toe. ak rize 
We can now return to the transformed hamiltonian and rewrite it taking into 


é Ay ^ ő ei? m A 
account (2.11), (2.14), (2.16a), (2.16b), (2.17b) and putting e^ - H, ,- e^ = H,,?). 
! ↄ pe. at E CiU menn 


4 A ^ n : 1 1 h 3. 
3) The problem of transforming en. H, dh ill be investigsted in paragrap 
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We get 
LS à seh c. EC NA coy So E E i = XVI „ 
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The expression (2.19c) describes a weak attractive electron — electron and 
ion — ion as well as repulsive electron — ion interactions. As Bohm and Pines have 
proved (1953b p. 620) this term may be neglected. In equation (2.19a) instead of m and 
M appear effective electronic and ion masses m* and M* given by the equality (2.19d). 
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Since the forces depending on polarisation of the ion — core are of short — 


. . a 1 3 A 
range because of their proportionality to — We can replace in H 
r T 


, the term of interac- 
tion among electrons and ions treated as points 


imo by more general expression including 
polarisation effects, i. e. 
3 de rij 
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In such a manner we obtain a hamiltonian describing better the electron behaviour 
in crystals, since V (|r; — Tal) represents now interaction energy including all the 


electrostatic forces acting among the electrons and ion — cores. 
J E 
3. The new subsidiary conditions 


Let us now investigate in what fashion the expression et'7'"i changes under 
the influence of canonical transformation. We have 
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and, in particular, 
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We get further 
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When we content ourselves with bai only two terms of expansion in series of 
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of this quantity with S (the same concerns the similar ion expression) we get 
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The terms AD, AD can be immediately obtained from the electron terms 
enclosed in the sam surly brackets by the substitution — e Ze, mM, n>N, 
r. . P >p. In the additional terms to H,, in the eqution (3.2c) the summing 
ranges over A and o do not coincide. Hence the additional terms are small. As to 
their physical meaning they cause the damping of the plasma oscillations. The 
neglect of them would mean that we do not take into account their disruptive 


influence on the organized behaviour of the plasma. We have then 
-iSh H asa. H 
ce Here (3.2d) 


In all calculations which led to (3.1f) we did not take into account that co; is an ope- 

rator which acts on er. Exact appreciation of this fact is almost impossible, because 
1 E 

we do not know the precise value of commutator „ z, e 
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Since we however restrict our calculations to first commutators with S, the error 


which arises from neglecting the operator character of w, is not large. 
We can proceed now to an exact analysis of the subsidiary conditions (2.3). 


For this purpose we write them in form 
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A suitable transformation cut off after the first commutators with S gives 
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including among others all higher range commutators. Making. use ae the random 
phase approximation and (3.1c), (3.1d) we get i 
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We apply to (3.4b) the dispersion relations (2.14) and finally obtain 
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When we omit the little quantity s (l, d, e , „„ D, D, 


and replace ein y by ꝙ we get new subsidiary conditions. The new state vector 
q loses its association with wave function Y and must be treated as a quantity to be 
determined anew from (3.4c). Since in the course of applying the canonical transfor- 
mation we have neglected small terms depending linearily on the operators a, and a, 
(they are hidden in e (a; , a}, .. ) and appear from higher, than first range, commuta- 
tors), it is not longer necessary to keep the term proportional to (1 — x,), especially 
as this term could be made as small as desired. For this purpose it is sufficient to 
take x, arbitrarily near to unity. Therefore we can write the subsidiary conditions 
in form 
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It seems to be quite sure that this level of Bohm- Pines theory must be treated 
rather as postulated one than correctly derived. 


When we assume that the ion — cores are approximately at rest and apply the 
adiabatic approximation we obtain the following final equations 
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I wish to thank to Professor S. Szczeniowski for valuable suggestions and helpful 


discussions during the writing of this article. 
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LETTERS TO THE EDITOR 


ON DESACCOMMODATION EFFECT IN NICKEL 
By Janusz MORKOWSKI 


Department of Ferromagnetics, Institute of Physics of the Polish Academy of Sciences 


(Received July 31, 1958) 


The author has observed a time decrease of reversible permeability after demagne- 
tization or stressing in specimens of Ni under compressive stress. Some preliminary 
results shall be presented here. 

In the former paper by the same author (sent to publication in Acta Physica 
Polonica) the apparatus for measuring reversible permeability under compressive 
stress was described. Experiments described here were performed on the same appara- 
tus. Momentary values of reversible permeability were determined from the values 
of e. m. f. induced in secondary coils of mutual inductances when the bridge was near 
balancing. The values of e. m. f. were measured directly on valve voltmeter with 
voltage amplifier at different moments of time or recorded on paper by means of 
an electronic recorder. 

The measurements were performed on long cylindrical specimens (0.48mm 
in diameter) from polycrystalline Ni (containing; Mg — 0.82%, Fe — 0.21%, Cu — 
0.08%, Si — 0.10%), annealed in hydrogen at about 1000°C for 2.5 hours and cooled 
to room temperature at a mean rate of about 12°C/minute. 

With regard to the smallness of the observed values of the time dependent part 
of the reversible permeability, the experimental errors in determining them are 
considerable (up to about 10%, of the value of u, (0) — u, (co), nevertheless some fa- 
cts may be established. 

The time decrease of permeability was observed merely in specimens at compres- 
sive stress. It is very probable that it exists also in another circumstances, but is 
then so small that it could not be detected in our experiments. 

After an increase of the value of compressive stress there could be observed 
a time decrease of the reversible permeability (measured at magnetic field strength 
0,45 oersted; all experiments were performed at this value of field) from the initial 
value, several percent ‘greater than its final stationary value — see Fig. 1. 
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Fig. 2 shows time decrease of the reversible permeability after demagnetization; 
differences between initial and stationary value of reversible permeability are of the 
order of 1% only. 

The effect observed depends strongly on temperature. On Fig. 3 there are pre- 
sented the results obtained after demagnetization at the same compressive stress but 
at various temperatures. 

The features of the described after effect show without doubt that it is identica 
with the so called „phenomenon of desaccommodation“ (see for instance Snoek 1947 
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F ig. 2. Time decrease of the reversible permeability after demagnetization at several values of compressive 
is stress at 16°C 


Becker, Döring 1939), observed in some ferromagnetic materials. A theoretical explana- 
tion of the desaccommodation effect was given by Snoek (1938, 1939) who supposes 
| that the desaccommodation is caused by gradual adaptation of 180? Bloch walls to 


itions i i ti mechanical shock. The process 
new positions in the course of time after magnetic or p 
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of adaptation is connected with the diffusion of C or N atoms, this should explain 
the strong temperature dependence of the desaccommodation. Néel (1951) pointed 
out some difficulties involved in Snoek’s theory and has shown theoretically that the 
desaccommodation is mainly connected with the adaptation of 90°-walls and that 
180°-walls can contribute to the desaccommodation in exceptional cases only. 

The assumption, however, that in the experiments described here the adapta- 
tion of 180°-walls which takes place leads to very simple qualitative explanation of 
the seemingly strange fact that desaccommodation is observed merely at compressive 
stress. It is well known that in Ni at compressive stresses the most important mechanism 
of magnetization changes consists in the displacements of 180°-walls. It is expected 
also that — in unstressed specimens — displacements of the 90°-walls determine 
the main part of the magnetization changes at small magnetic field strength overriding 
the influence of 180°-walls (cf. Vonsovski, Shur 1948, p. 359). Thus the time dependent 
part of the reversible permeability — determined as we have assumed by-adaptation 
of 180°-walls — is too small, as compared with the stationary one, to be detected 
in our experiments with unstressed specimens, being significant only when the 180°- 
walls displacements prevail. Following similar arguments the absence of the desaccom- 
modation effect at tensile stresses may be explained. 

I wish to thank Prof. Dr. Sz. Szczeniowski for his interest in this work. 
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APPLICATION OF THE FINE FOCUS X-RAY TUBE FOR THE 
INVESTIGATION OF THE SIZE AND ORIENTATION AND RELATIVE 
POSITION OF THE BLOCKS OF MOSAIC IN MONOCRYSTALS 


By J. AULEYTNER 


(Received August 16 1958) 


In this paper are presented further possibilities of a method described previously 
by the author (Auleytner 1958), in particular, the simultaneous determination of the size 
of blocks of the mosaic and their twists (50 u, 30°). It is possible in certain cases to 
determine the size of blocks of the order of several microns. The method permits the 
analysis of the relative position of the blocks, which can be of great importance in 
investigations of the dislocation structure of crystals. 


In a previous paper by the author (Auleytner 1958) there was given a method for 
determining the disorientation of the blocks of the mosaic, which made possible 
a detailed investigation of crystal surfaces. Discussed below are further possibilities 
of this method, in particular, the determination of the size of regions of different orien- 
tation and the investigation of their positions. 

As already mentioned in a previous paper, if the crystal is displaced in the XY 
or YX directions (Fig. 1) by a segment greater or equal to the width of the illuminated 


region one may obtain, in a relatively short time, a picture of the mosaic structure of 


a larger area of the monocrystal surface. On a single photographic plate one obtains 
a number of diffraction lines, which gives an image of the structure of the individual 
parts of the surface of the crystal being investigated in detail. Fig 2 may serve as an 
example. In it are shown the diffraction lines K, K,, Cu obtained by the above method 
in the case of a zinc monocrystal (the investigated cleavage surface of the crystal coinci- 
- ded with the diffraction plane (0001); the width of the illuminated region was fixed at 
| ; (81) 


_ 
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0.05 cm., the magnification x 4). Considering the shape of the equivalent parts of the 
lines on the level denoted by an arrow, we note that after shifting the crystal by 0.2 mm 
(lines 4, 5, 6, 7), the shape of line (8) changes. This indicates the passage to a region of dif- 
ferent orientation whose width is 0.2 mm. An additional shift by 0.1 mm reveals further 
passage to a region of still another orientation. Analysing the diffraction lines on the 
level of arrow B in Fig 3, we note that the twisting of one block (line 1—2) is related 
to the twisting of the following blocks (lines 3—4, 5—6, 7). A similar coupling of the 
twist of blocks was observed in the case of a cleavage surface of K Br crystals. Fig. 4 
shows a comparison of the K, K, Cu lines obtained for zinc monocrystals in the case 
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2) linear dimensions of the blocks smaller than the width of the illuminated region. 

In the first case, when the beam moves over the area of one block the width of the 
diffraction line and the angular position remains constant. At the boundary of two 
blocks of different orientation, when the incident beam begins to move into the next 
block there appears either an additional diffraction line, or the width of the diffraction 
line increases, becoming asymmetric (Fig. 3, line 3 on the level of arrow 4). The resol- 
ving power in this case depends on the accuracy of setting the surface of the crystal 
on the axis of rotation of the spectrometer table, the divergence of the incident beam, 
the characteristic line’s width for a given radiation and the type of crystal (Bac- 
kovsky 1956, 1957). The measurement of the angle of twist is made by means of 
a comparator where as a point of reference is taken a reflex whose shape or position 
which does not change in passing from one line to another. The accuracy of the measur- 
ement is determined by the control in shifting the crystal. 

In the second case the diffraction line can consist of many reflexes displace with 
respect to one another. If the angles of twist of the blocks are smaller than the angular 
widths of the individual reflexes, then the diffraction line is diffuse. The diffuseness 
appears also as a result of the large dispersion of the size of the blocks. If the size of the 
blocks and their mutual twist were uniform, then, in principle, it would be possible 
to observe the separation of the diffraction lines coming from blocks of linear dimen- 
sions of the order of several microns. 

Conclusions: The above method permits the simultaneous determination of the 
size of blocks of the mosaic and their twists of the order of 504 and 30”. It is possible 
in certain cases to determine the size of blocks of the order of several microns. This 
method permits the analysis of the relative position of the blocks, which can be of 
great importance in investigations of the dislocation structure and plastic deformation 
of crystals. 

In employing this method for investigating the mosaic structure of Zn and KBr 
crystals, it was noted that the mutual twisting of two blocks is connected with the 
twisting of a number of neighbouring blocks. 

I would like to express my sincere gratitude to Professor Dr L. Sosnowski for 
his valuable remarks and constant interest. 
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